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1.  INTRODUCTION 

Work  in  our  laboratory  and  others  shows  that  the  energy  landscape  of  a  protein  is 
complex  consisting  of  a  large  number  of  energy  valleys  separated  by  energy  mountains  and 
ridges.  Moreover,  experimental  and  theoretical  evidence  is  accumulating  that  the  energy 
landscape  is  hierarchical  and  contains  valleys  within  valleys.  We  continue  our  investigation 
of  several  aspects  of  protein  dynamics;  (i)  the  structure  and  organization  of  the  energy 
landscape  focusing  on  the  hierarchical  nature  of  the  landscape,  (ii)  protein  motions  within 
the  energy  landscape,  (iii)  the  laws  governing  the  motions,  and  (iv)  the  slaved  glass 
transition  in  proteins.  We  also  attempt  to  uncover  the  laws  of  protein  reactions  under 
pressure  (reaction  theory). 

We  study  phenomena  in  detail  in  one  protein  (myoglobin)  and  then  explore  how  the 
phenomena  depend  on  protein  structure.  To  do  this  we  move  to  other,  possibly  more 
complicated,  proteins  and  to  synthetic  mutant  myoglobins.  Our  experimental  approach 
involves  flash  photolysis  applied  over  wide  ranges  in  temperatures  (10  to  330K),  time 
(ns  to  ks),  and  pressure  (0.1  to  300  MPa).  Pressure  is  central  to  our  work  both  as  static 
parameter  and  dynamic  perturbation  to  probe  protein  relaxations.  We  monitor  protein 
reactions  and  relaxations  from  the  near  ultraviolet  to  the  mid-infrared  with  Fourier 
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Transform  Infrared  Spectroscopy  playing  an  essential  role.  We  also  use  theoretical  and 
computational  tools  to  build  models  of  protein  dynamics  and  reactions. 

2.  ACCOMPLISHMENTS 


Progress  has  been  very  good  during  the  first  year.  A  short  paper  entitled  "Glassy 
Behavior  of  a  Protein"  has  been  published  which  shows  that  proteins  and  glasses  share 
essential  characteristics  (1).  In  particular  we  proved  that  myoglobin,  like  glasses,  exhibits 
metastability  below  a  transition  temperature  and  relaxation  processes  that  are 
nonexponential  in  time  and  non- Arrhenius  in  temperature.  A  second,  much  longer, 
manuscript  entided  "Proteins  and  Pressure"  has  appeared  in  the  Journal  of  Physical 
Chemistry  (2).  This  manuscript  describes  many  of  our  discoveries  regarding  protein 
relaxations  and  reactions  obtained  using  wide  ranges  of  pressure  and  temperature.  These 
results  greatly  clarify  the  effects  of  pressure  and  imply  that  most  previous  work  on  pressure 
effects  in  proteins  needs  to  be  reexamined.  Both  papers  are  included  in  the  Appendix  to 
this  Report.  Our  work  also  includes  the  following  activities: 


(i)  Pressure  Effects  on  CO  Binding  to  Myoglobin.  We  performed  a  series  of 
experiments  using  FOR  spectroscopy  to  probe  the  effects  of  pressure  on  the  kinetics  of  CO 
binding  to  the  different  A  substates  of  myoglobin  after  photodissociation  below  lOOK.  The 
A  substates  involve  tier  0  of  the  hierarchy  of  substates.  We  use  different  pathways  in  the 
pressure-temperature  plane  to  separate  conformational  (equilibrium)  and  kinetic 
phenomena.  These  experiments  clearly  demonstrate  the  subtle  and  important  effects  of 


pressure  on  protein  dynamics  and  function.  Our  main  rebu’fs '  ’pear  in  the  paper  "Proteins 
and  Pressure"  (2). 


in  For 


(ii)  Pressure-Temperature  Studies  of  Horse  and  Sperm  Whale  Myoglobin.  The  ,0,3 

kinetics  of  horse  and  sperm  whale  myoglobin  are  very  similar  at  atmospheric  pressure  even  . 

when  temperature  is  varied  over  a  wide  range.  However  studies  on  tiitae  two  mamimlian - 

■-  J  on/ 

myoglobins  over  wide  ranges  of  temperature  and  pressure  reveal  significant  differences  i  Codes 

^  ..ii  and/or 

.-'-lat  I  ''paclal 


□  □ 
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which  are  mainly  a  result  of  conformational  effects  (2).  It  is  interesting  to  speculate  if  these 
differences  are  related  to  the  fact  Uiat  whales  must  dive  but  horses  do  not! 

(iii)  Quasistatic  and  Kinetic  Studies  of  More  Complicated  Heme  Proteins.  VVe  have 
begun  to  study  the  effects  of  protein  structure  using  the  more  complicated  heme  proteins 
horseradish  peroxidase  type  C  (HRP-C)  and  cytochrome  P450.  We  use  CO  binding  to  the 
heme  iron  in  HRP-C  and  P450  to  probe  the  slaved  glass  transition  and  protein  relaxation 
processes.  Using  pressure,  we  have  observed  CO  stretch  bands  in  HRP-C-CO  which 
correspond  to  protein  structures  not  previously  resolved.  P450  with  substrate  bound  and 
unbound  has  also  yielded  surprising  results.  We  have  begun  to  study  the  complex 
spectrum  of  relaxation  processes  in  these  proteins  near  the  slaved  glass  transition  and  have 
presented  preliminary  reports  of  our  results  at  the  March  Meeting  of  the  American  Physical 
Society  (3)  and  the  Meeting  of  the  Biophysical  Society  (4).  We  will  continue  our  studies 
concentrating  on  samples  and  conditions  of  temperature  and  pH  where  the  pressure  effects 
are  particularly  dramatic.  This  work  may  give  significant  new  information  regarding  the 
substates  of  tier  0  (CS®)  which  we  believe  are  crucial  to  protein  function. 

(iv)  General  Aspects  -  Pressure  and  Proteins.  Our  previous  work  has  shown  that  the 
substates  of  tier  0  in  myoglobin  (A  substates)  have  different  kinetics  for  rebinding  CO  and 
are  affected  differently  by  pressure.  The  manuscript  "Proteins  and  Pressure"  (2)  also 
includes  a  reexamination  of  our  previous  work,  much  of  it  unpublished,  on  pressure  effects 
in  proteins.  We  analyze  in  detail  the  subtle  but  crucial  effects  of  the  different  protein 
structures  corresponding  to  tier  0  (CS®)  of  the  hierarchy  of  substates  in  myoglobin.  Since 
the  different  substates  of  CS^  also  bind  CO  at  different  rates  and  have  different  structures 
and  volumes,  the  switch  from  one  substate  to  another  provides  the  possibility  of  an 
efficient  control  mechanism.  We  have  shown  in  Ref.  (2)  that  the  effects  of  pressure  in  the 
visible  spectral  region  can  be  very  misleading  since  this  spectral  region  is  not  sufficiently 
sensiL’.  c  to  conformauonal  effects. 


4 


3.  CONCLUSIONS 

The  technique  and  results  under  this  contract  may  be  important  for  both  biology  and 
physics.  For  biology,  detailed  knowledge  of  protein  motions  at  different  time  and  length 
scales  is  necessary  for  an  understanding  of  protein  and  enzyme  reactions  at  the  molecular 
level.  In  particular  the  different  behavior  between  horse  and  whale  myoglobin  under 
pressure  is  intriguing  and  calls  for  further  studies.  For  physics,  proteins  may  well  become 
paradigms  of  complex  systems  including  glasses,  spin  glasses,  and  neural  networks. 
References  (1)  and  (2)  can  be  consulted  for  additional  details  and  discussion. 
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Quasisialic  and  kinetic  studies  of  the  infrared  CO  stretch  bands  of  carbonmonoxymyoglobin  show  that 
proteins  and  glasses  share  essential  characteristics,  in  particular  metastability  below  a  transition  temper¬ 
ature  and  relaxation  processes  that  are  nonexponential  in  time  and  non-Arrhenius  in  temperature. 


PACS  numbers  87.15  Da,  64.70  Pf,  78.30.Jw.,  82  20.Rp 

Proteins  and  glasses  may  appear  to  have  little  in  com¬ 
mon.  Proteins  are  macromolecules  with  well  defined 
structures';  glasses  are  frozen  liquids.  Despite  this 
difference,  proteins  and  glasses  share  one  fundamental 
property:  Both  can  assume  a  very  large  number  of  near¬ 
ly  isoenergetic  conformational  substates  (CS),  valleys  in 
the  conformational  energy  landscape.  For  proteins,  the 
existence  of  CS  followed  from  the  nonexponential  time 
dependence  of  the  binding  of  small  molecules  (Oi  and 
CO)  to  myoglobin  at  low  temperatures.^  Supporting  evi¬ 
dence  came  from  other  experiments’  and  from  theory.'* 
For  glasses,  a  potential-energy  surface  with  a  large  num¬ 
ber  of  minima  was  postulated  by  Goldstein’;  for  spin 
glasses,  the  evidence  came  from  theory.*  The  existence 
of  CS  in  proteins  and  glasses  raises  the  question  as  to 
whether  these  systems  share  other  properties.  We  now 
describe  some  attributes  of  glasses  and  later  show  that 
these  are  also  found  in  proteins.^ 

Glasses  are  formed  when,  on  cooling,  a  liquid  becomes 
a  structurally  disordered  solid.*  The  temperature  at 
which  the  viscosity  reaches  10'’  poise  is  called  the  glass 
temperature,  Tg.  The  specific  heat  below  1  K  is  approxi¬ 
mately  proportional  to  the  temperature.’  Glass  proper¬ 
ties  well  below  Tg  depend  on  history;  glasses  are  in  a 
metastable  (nonequilibrium)  state.  Near  and  above  Tg 
the  response  of  a  glass  to  a  mechanical  or  electrical  per¬ 
turbation  is  dominated  by  the  a  relaxation.  Its  relaxa¬ 
tion  function  <l>r(f)  is  usually  nonexponential  in  time  and 
can  be  parametrized  by  a  stretched  exponential,  <t>rU) 

— exp[-  (it/)^!,  or  by  a  power  law, 

d>,(/)-fl+K.(r)fI"".  (1) 

The  average  rate  at  temperature  T  is  (k)  “nkriT).  The 
temperature  dependence  of  kriT)  follows  the  Arrhenius 
relation,  /t,(r)  “/4exp[  — £/A:fl7'l,  only  over  small  tem¬ 
perature  intervals.  Typical  values  near  Tg,  £  ~  1.6  eV, 
/I“I0'"  s~',  also  imply  that  the  Arrhenius  relation  is 
inappropriate  for  glasses.  However,  kr(T)  can  be  de¬ 
scribed  over  more  than  10  orders  of  magnitude  either  by 
the  Vogel-Tammann-Fulcher  equation,* 

it.(7-)-/lvTFCxp[-£/i(fl(r-ro)l  (2) 


or  by  the  relation 

kr(T)-koexpl-(To/T)"].  (3) 

Both  relations  fit  the  data  for  glycerol  (TpaslSS  K) 
from  190  to  260  K.''’-'^ 

We  now  examine  the  proteins  for  glasslike  properties. 
Two  are  well  known;  Each  individual  protein  is  disor¬ 
dered  (aperiodic)  and  the  specific  heat  of  proteins  below 
1  K  is  glasslike.'’  The  other  attributes,  however,  have 
been  less  well  explored.  Here  we  report  experiments  that 
verify  the  metastability  at  low  temperatures,  and  the 
nonexponential  time  and  the  non-Arrhenius  temperature 
dependence  of  the  protein  relaxations  near  200  K  in  car¬ 
bonmonoxymyoglobin  (MbCO). 

The  folded  polypeptide  chain  of  the  oxygen-storage 
protein  myoglobin  (Mb)  embeds  a  heme  group  with  a 
central  iron  atom  which  reversibly  binds  ligands  such  as 
O’  and  CO. '  Our  experiments  focus  on  the  stretch 
bands  of  CO  bound  to  Mb  which  are  very  sensitive  to 
external  parameters  such  as  solvent,  pH,  temperature 
(r),  and  pressure  (£).'*■"  We  measure  the  stretch 
bands  with  a  Mattson  Fourier  transform  infrareo  spec¬ 
trometer.  Figure  1  shows  that  MbCO  displays  at  least 
three  different  CO  stretch  bands,  Ao,  At,  and  A}.  Fits  to 
Voigtian  line  shapes'*  yield  the  areas  (A,),  center  fre¬ 
quencies  (v,),  and  linewidths  (F,)  of  the  A  bands.  We 
also  observe  the  rate  of  heat  absorption  via  differential 
scanning  calorimetry  (DSC).  The  experiments  fall  into 
two  classes,  quasistatic  and  kinetic.  Quasistatic  indicates 
that  the  glasslike  behavior  of  MbCO  below  a  transition 
temperature  Tsg  prevents  attainment  of  thermodynamic 
equilibrium.  Quasistatic  measurements  determine  the 
band  parameters  as  functions  of  solvent,  pH,  T,  and  P. 
In  the  kinetic  studies  we  observe  the  relaxation  of  the 
protein  after  a  pressure  release. 

Quasisialic  experimenis. —  Figure  1  shows  the  ir  spjec- 
tra  from  1910  to  1990  cm  “'  in  a  75%  glycerol-water  sol¬ 
vent  (3:1  by  volume)  at  pH  6.8  with  potassium  phos¬ 
phate  buffer.  The  sample  was  brought  to  a  pressure  P  at 
300  K.  Data  were  then  taken  under  constant  pressure  at 
successively  lower  temperatures.  The  cooling  rate  of 
0.01  K/s  and  waiting  time  of  about  600  s  at  each  temper- 
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FIG.  I.  Absorption  spectra  of  the  CO  stretch  bands  of 
MbCO  at  100  MPa,  Solvent:  75%  glycerol  and  water  mixture, 
buffered  with  100  mM  potassium  phosphate  to  pH  6.8  Inset: 
Cross  section  through  the  heme  showing  the  angle  a,  beiween 
the  CO  dipole  and  the  heme  normal. 

ature  were  kept  constant.  Measurement  of  an  ir  spec¬ 
trum  took  about  250  s.  Figure  2  gives  the  temperature 
dependence  of  the  ratio  ro^A{)lA\  at  three  different 
pressures.  The  ratio  /-q  exhibits  five  regions:  (i)  Glass 
region.  The  protein  is  in  a  metastable  state  in  which 
large-scale  motions  are  too  slew  to  be  observed.  We 
prove  metastabilily  by  showing  that  the  ir  spectrum  de¬ 
pends  strongly  on  the  path  by  which  a  point  in  the  (T,P) 
plane  is  reached,  (ii)  Slated  glass  transition.  The  bend 
at  Esg  in  Fig.  2  reveals  the  transition  from  an  equilibri¬ 
um  situation  well  above  Fsg  to  a  metastable  one  well 
below.  A  similar  bend  occurs  in  r(7'),'''  and  the  transi¬ 
tion  is  also  visible  in  DSC.  The  glass  temperature  T„g  of 
the  protein  depends  crucially  on  the  glass  temperature 
Tg  of  the  solvent;  so  we  call  the  transition  “slaved”''*: 
DSC  data  at  a  heating  rate  of  0.083  K/s  in  75%  glycerol 
and  water  mixture  give  rg“175  K  and  Tsg”178  K;  in 
60%  ethylene-glycol-water  mixture  Fg”l49  K  and 
F,g“*155  K;  in  solid  poly(vinyl)alcohol  Fg  >  280  K  and 
Tig>  280  K.  Fsg  and  Tg  depend  approximately  logarith¬ 
mically  on  cooling  or  heating  rate;  a  factor  of  10  increase 
in  heating  rate  raises  Fsg  by  about  5  K.  Fsg  obtained 
from  the  bend  in  /-q  is  about  20  K  higher  than  that  from 
DSC.  With  increasing  pressure  Fsg  increases  weakly  as 
is  apparent  from  Fig.  2.  (iii)  Lower  equilibrium  region. 
Transitions  among  the  A  substates  occur  faster  than  the 
time  of  observation  so  that  the  A  substates  are  in 
thermal  equilibrium.  The  temperature  and  pressure 
dependences  of  the  ratios  r,  determine  the  relative  ener¬ 
gies,  entropies,  and  volumes  of  the  A  substates.  (iv) 
Upper  equilibrium  region.  From  250-320  K,  the  ratio 
ro  shows  a  remarkable  reversal,  (v)  Unfolding.  The 
protein  begins  to  unfold  irreversibly.  Since  regions  (iv) 
and  (v)  are  not  the  focus  of  this  Letter  we  give  no  fur- 
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FIG.  2.  Plot  o(  ro~.4olA\  vs  IOVF  for  three  pressures  Ro¬ 
man  numerals  signify  the  five  temperature  regions  defined  in 
the  text.  The  transition  temperature  F,g  is  indicated.  The 
dashed  line  shows  the  extrapolation  from  region  (lii)  for  1 
MPa. 

ther  details. 

Relaxation  experiments. —  In  a  pressure-release  exper¬ 
iment,  100  MPa  is  applied  at  240  K,  the  sample  is  cooled 
to  the  desired  temperature  and  held  there  for  600  s  fol¬ 
lowed  by  measurement  of  an  ir  spectrum  for  about  250  s. 
The  pressure  is  then  released  to  7  MPa  in  a  few  seconds 
and  ir  spectra  are  taken  at  approximately  exponentially 
increasing  times  from  about  10  to  more  than  10*  s.  The 
time-dependent  behavior  of  the  A  bands  after  pressure 
release  reveals  several  relaxation  processes.  To  interpret 
the  relaxation  data,  we  assume  that  the  conformational 
substates  of  a  protein  are  arranged  in  a  hierarchy  of  at 
least  two  tiers,  CS"  and  CS'.'*"*  The  different  stretch 
bands  correspond  to  substates  of  tier  0  (CS°).  We 
characterize  each  substate  A,  by  the  center  frequency  v', 
of  the  CO  stretch  band  (Fig.  1)  and  by  the  angle  a,  be¬ 
tween  the  CO  dipole  and  the  heme  normal  (inset  of  Fig. 
1).'’  It  is  likely  that  the  substates  A,  differ  not  only  in 
the  angle  a,,  but  also  in  the  overall  structure  of  the  mole¬ 
cule.  A  protein  may  move  from  one  substate  of  tier  0, 
say  Ao.  to  another,  say  A\.  We  call  this  motion  FIM  0, 
where  FIM  is  an  acronym  for  “functionally  important 
motion,”'"  Rebinding  of  CO  after  flash  photolysis  to 
each  of  the  ,4  substates  is  nonexponential  in  time.'*  This 
observation  implies  that  each  substate  of  tier  0  is  subdi¬ 
vided  into  substates  of  tier  I,  CS'.  We  denote  as  FIM  1 
the  nonequilibrium  conformational  motions  between  sub¬ 
states  of  tier  I.  With  these  concepts  we  explain  our  ex¬ 
periments. 

Below  160  R — The  center  frequencies  of  the  A  bands 
shift  (0.4  cm~'  for  Ao)  faster  than  our  shortest  mea¬ 
surement  time  (10  s),  but  no  changes  in  relative  intensi¬ 
ties  occur  up  to  10*  s.  Since  the  fast  shift  occurs  even  at 
15  K,  we  interpret  it  as  an  elastic  relaxation. 

165-IS5  K.  —  The  area  of  Ao  remains  constant  to 
within  2%  from  10  to  more  than  10*  s  while  the  peak 
narrows  and  its  amplitude  increases.  The  center  fre- 
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quencv  vo  ir  y  shifts  rapidly  from  1964.5  to  1964.9 
cm  ~  ‘  and  then  moves  measurably  towards  its  low- 
pressure  value,  1966.1  cm  The  shift  in  vq  can  be  de¬ 
scribed  by  a  relaxation  function  <l>i (t)  —  [v(|(/)  —  v',  (oo)]/ 
[vo(0'^  1  —  V()(oo)],  where  i  is  the  time  after  pressure 
release  and  vofO'*’)  represents  the  center  frequency  im¬ 
mediately  after  the  rapid  elastic  relaxation.  For 
we  use  two  approximations,  the  7  MPa  value  from  a  qua¬ 
sistatic  experiment  at  the  same  temperature  and  the 
value  obtained  by  extrapolation  from  region  (iii).  Figure 
3  shows  the  relaxation  function  Olf/)  using  the  former 
choice  for  vofoo).  'I>|(/)  is  not  greatly  changed  with  the 
other  choice.  We  interpret  the  relaxation  characterized 
by  <t>i  as  FIM  1,  a  redistribution  of  substates  CS'  within 
the  substate  Ao- 

To  evaluate  the  data  in  Fig.  3,  we  parametrize  <hi(r) 
by  the  power-law  equation  (1).'°  If  we  further  assume 
that  kr(T)  is  determined  by  the  Arrhenius  law,  we  ob- 
tain  £  “  1 .3  eV, /I  ”  10^'’ s  These  values  are  similar 
to  the  values  quoted  above  for  the  a  relaxation  in  glasses 
near  they  imply  that  the  Arrhenius  relation  is  inap¬ 
propriate  for  describing  FIM  1  and  suggest  the  use  of 
Eq.  (2)  or  (3).  We  select  Eq.  (3)  to  fit  kr(T)  because  it 
uses  only  two  parameters  and  therefore  provides  less  am¬ 
biguous  extrapolations.  A  fit  with  Eqs.  (I)  and  (3),  us¬ 
ing  a  Levenberg-Marquardt  algorithm  to  minimize  x'<  is 
shown  as  solid  lines  in  Fig.  3  and  yields  log(<:o/s“') 
-1 7.0 ±2.5,  To- 1130 ±80  K,  n -0.26  ± 0. 1 2.  The 
errors  are  conservative.  We  have  also  measured  the  a 
relaxation  in  75%  glycerol  and  water  mixture  using 
specific-heat  spectroscopy and  obtain  log(An/s  ” ' )  «  18 
and  ro«  1130  K.  The  values  of  ko  and  To  for  the  pro¬ 
tein  and  the  solvent  are  remarkably  close  and  support  the 
notion  of  a  slaved-glass  transition:  Not  only  do  the  .sol¬ 
vent  and  protein  exhibit  comparable  glass  temperatures, 
the  a  relaxation  of  the  solvent  and  FIM  1  in  the  protein 
obey  nearly  identical  temperature  dependences. 

A  second  relaxation,  no^ shown  here,  occurs  between 
170  and  190  K:  A\  and  Ay  interconvert,  with  A\  de¬ 
creasing  and  Ay  increasing  nonexponentially  in  time;  the 


FIG.  3.  Relaxation  function  <t>i(r)  for  the  center  frequency 
Vo  of  the  Aa  band  vs  time.  Solid  lines  are  fits  with  Eqs.  ( 1 )  and 
(3)  with  parameters  as  given  in  the  text. 


barrier  between  A\  and  Ay  is  smaller  than  that  between 
Af)  and  the  other  A  substates."' 

190-210  K. —  FIM  1  and  the  interconversion  between 
A]  and  Ay  are  faster  than  10  s  so  that  after  pressure 
release  the  center  frequencies,  the  linewidths,  and  the  ra¬ 
tio  Ay/A  \  have  the  values  measured  quasistatically  at  7 
MPa.  A  new  relaxation  process  is  observed,  the  ex¬ 
change  of  Ao  with  A]+Ay.  We  interpret  this  exchange 
among  substates  of  tier  0  as  a  conformational  change  of 
the  entire  protein  and  call  it  FIM  0.  We  define  a  re¬ 
laxation  function  for  FIM  0  through  <I>o(/ )“[/!()(/ ) 
— /I(i(<»)]/[zlf)(0''^ )  — The  area  is 

found  by  extrapolation  from  region  (iii)  as  indicated  by 
the  dashed  line  in  Fig.  2."'  <t>o(t),  shown  in  Fig.  4, 
displays  more  structure  than  FIM  1  and  a  simple  power 
law  or  stretched  exponential  does  not  fit  the  data  well. 
Nevertheless,  insight  is  gained  by  fitting  the  data  with 
Eqs.  (1)  and  (3),  with  the  result  log(A.o/s  ~ ~  1 4,  7o 
=  1200  K,  n  5=0.8.  The  fit  reproduces  the  general  be¬ 
havior  of  but  the  detailed  structure  will  require  a 

more  elaborate  model. 

Our  results  permit  four  conclusions:  (1)  Proteins  and 
glasses  indeed  share  common  properties  in  addition  to 
the  existence  of  many  valleys  (conformational  substates) 
in  the  energy  landscape.  The  similarities  include  meta¬ 
stability  below  the  glass  temperature  and  the  presence  of 
a  relaxation  process  governing  large-scale  motions  (a  re¬ 
laxation  in  glasses,  FIM  1  in  proteins)  which  is  nonex¬ 
ponential  in  time  and  does  not  follow  the  Arrhenius  law, 
(2)  While  glasses  usually  show  only  one  relaxation  pro¬ 
cess  with  the  general  characteristics  of  the  a  relaxation, 
MbCO  shows  at  least  two,  FIM  0  and  FIM  I.  FIM  0  is 
slower  than  FIM  I  and  probably  describes  overall 
changes  in  the  entire  protein,  cha.-acterized  by  different 
CO  stretch  frequencies  and  angles  between  the  bound 
CO  and  the  heme  normal.  FIM  1  describes  large-scale 
motions  within  a  given  substate  CS”.  (3)  The  protein  re¬ 
laxation  process  FIM  1  is  slaved  to  the  solvent;  the  tran¬ 
sition  temperature  Tsg  in  the  protein  is  close  to  the  glass 
transition  of  the  solvent.  Moreover,  FIM  1  in  the  protein 


FIG.  4.  Relaxation  function  <l>o(()  for  the  area  Ao  vs  time. 
Solid  lines  are  fits  with  Eqs.  (I)  and  (3)  with  parameters  as 
given  m  the  text. 
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and  the  a  relaxation  in  the  solvent,  at  least  in  the  case  of 
glycerol  and  water,  possess  remarkably  similar  tempera¬ 
ture  and  time  dependences.  The  hydration  shell  may 
well  play  an  important  role  in  this  coupling  between  pro¬ 
tein  and  solvent.*’  (4)  FIM  0  and  FIM  1  can  be  extra¬ 
polated  to  300  K  with  Eq.  (3)  with  the  result  <A;)=  10' 
s“'  for  FIM  0  and  (A:>~10"’s‘'  for  FIM  1.  While 
these  extrapolations  are  speculative,  they  indicate  where 
to  look  experimentally  and  suggest  that  molecular  dy¬ 
namics’’  may  be  able  to  simulate  FIM  1  hut  probably 
cannot  yet  reach  FIM  0. 

The  technique  and  the  results  presented  here  may  be 
important  for  both  physics  and  biology.  For  physics, 
proteins  may  well  become  paradigms  of  complex  sys¬ 
tems.  The  combination  of  the  pressure  jump  appro.ich 
with  site-specihe  spectroscopic  observation  permits  site- 
spccitic  studies  of  relaxation  phenomena  in  amorphous 
systems  When  applied  to  genetically  engineered  pro¬ 
teins,’  the  pressure  jump  technique  would  probe  the  re¬ 
laxation  of  specifically  designed  systems.  For  biology, 
detailed  knowledge  of  protein  motions  at  dilferent  time 
and  length  scales  is  necessary  for  an  understanding  of 
protein  and  enzyme  reactions  at  the  molecular  level. 
Slaving  may  be  an  efficient  control  mechanism  in  cells 
and  membranes. 
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The  exploration  of  the  combined  effects  of  pressure  and  temperature  leads  to  new  insights  into  the  structure,  dynamics,  and 
reactions  of  proteins.  A  protein  molecule  px)sses,ses  many  conformational  substates,  nearly  isoenergetic  structures  performing 
the  same  function  but  possibly  with  different  rates  Pressure  a'iers  both  the  relative  substatc  populations  and  the  functional 
properties  of  the  individual  substates  and  thus  can  produce  large  changes  in  the  static  and  dynamic  properties  of  a  protein 
.•nscmble  Such  prcssuie  effects  have  been  measured  for  myoglobin  (.Mb),  using  visible  and  FTIR  techniques.  .Above  a 
glass  temperature,  r,j,  pressure  redistriButes  substate  populations  and  c,.anges  protein  structure,  the  position  of  spectral 
bands,  and  protein  reaction  rates.  E-iquilbrium  substate  populations  measured  from  230  to  330  K  at  pressures  up  to  200  MPa 
vield  a  thermodynamic  comparison  of  different  substates  at  pH  5  '■  and  6  6  Near  T,^.  pressure  ju  p  experiments  reveal 
■.everal  relaxation  pnxesses.  Well  below'  the  properties  of  the  proteins  at  temperature  T  and  pressuie  P  depend  on  history. 
1  e  ,  on  the  order  in  which  the  sample  is  pressurized  and  cooled  The  effect  of  pressure  on  the  bi.iding  CO  and  Oj  to  Mb 
after  photodissociation  points  to  the  strong  influence  of  conformational  substates  on  protein  reactions.  Pressure  controls 
the  rate  of  binding  through  the  activation  volume  and  through  redistribution  of  substate  populations.  The  conformational 
effect  of  pressure  differs  for  CO  and  Oj.  and  for  sperm  whale  and  horse  mycglobin 


I  Introduction 

The  first  draft  of  this  paper  was  written  1 1  years  ago.  W'e  did 
not  publish  then  because  we  did  not  sufficiently  understand  the 
result.s,  but  we  continued  tc  explore  the  effects  of  pressure  on 
protein  slates  and  reactions.  During  this  time.  Professor  Harry 
Drickamer  has  been  an  unfailing  advisor  and  fnend  who  has  helped 
us  with  the  design  of  equipment,  with  overcoming  many  difficult 
problems,  and  most  importantly  by  setting  standards  for  scientific 
research.  We  now  believe  that  we  have  come  closer  to  compre¬ 
hending  the  effect  of  pressure  on  proteins,  and  we  u.sc  the  occasion 
of  the  Drickamer  Special  Issue  to  summarize  our  work.  Many 
collaborators  have  contributed  during  these  past  1 1  years,  ex¬ 
plaining  the  length  of  the  author  list  which  would  be  more  ap¬ 
propriate  for  a  high-energy  physics  paper  The  many  years  of 
work  are  responsible  fo'  the  length  of  the  paper. 

Pressure  is  an  important  variable  not  only  in  chemistry  and 
physics,  but  also  in  biology.'  Pressure  can  deactivate  enzymes, 
viruses,  and  toxins,  kill  bacteria  and  yeasts,  reverse  narcosis  in¬ 
duced  by  alcohol,  and  influence  bioluminescence.’  Sharks 
function  from  sea  level  to  extreme  depth,’-*  encountering  pressures 
from  0  1  to  100  MPa  ’  In  ord"r  to  explain  these  biological 
phenomena  on  a  molecular  'evel,  the  effects  of  pressure  on  isolated 
proteins  must  be  understood.  Many  pressure  studies  of  proteins 
have  been  performed,'  but  usually  in  a  narrow  tempe-ature 
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range  near  300  K  As  we  shall  show,  studies  over  a  wide  tem¬ 
perature  range  are  needed  for  an  interpretation  of  the  observed 
pressure  effects 

The  effects  of  a  pressure  P  on  a  protein  can  be  explained  with 
Figure  1  which  gives  t.vo  vastly  simplified  cross  sections  through 
the  energy  hypersurface  of  a  protein.  We  assume  th  '  the  protein 
has  two  conformations  (or  conformational  substates)  0  and  1 
t Figure  la)  and  two  states  .'X  and  B  (Figure  lb).  In  each  of  the 
two  substates  0  and  1 ,  the  protein  can  make  the  transition  B 
A.  Figure  I  shows  in  panel  a  the  cross  section  through  the  energy 
surface  at  constant  reaction  coordinate,  say  rc  =  B.  Panel  b 
represents  the  cross  section  at  constant  conformational  coordinate 
(say  cc  =  1 ).  Pressure  affects  the  protein  characterized  by  Figure 
1  in  two  ways.  In  equilihnum.  the  ratio  "g/w,  of  the  populations 
H'o  and  H  ,  of  the  two  conformational  substates  0  and  1  at  tem¬ 
perature  T  is  given  by 

hq/w,  =  exp[-Af7//?r]  (1) 

where  AC  is  the  difference  in  Gibbs  free  energy  between  the  two 
substates. 

AG  =  A£ -F  P  AF- TAA  (2) 

Here  Af.  AF.  and  A5  are  the  differences  in  internal  energy, 
volume,  and  entropy  between  the  two  conformations.  Pressure 
shifts  the  equilibrium  ratio  by  the  factor  exp[-P  AF].  Pressure 
also  affects  the  reaction  rate.  If  a  reaction  P  -*  A  occurs  in  the 
conformational  substate  i,  the  reaction  rate  coefficient  can 
F  written  as 

=  e  exp(-G*B^/RT)  (3) 

Here 

G*BA  =  f'BA  +  /’^’•BA  -  TS\^  (4) 

is  the  activation  Gibbs  free  energy  and  v  is  the  frequency  factor 
in  the  substate  i.  The  star  indicates  that  G*  is  ideally  measured 
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Figure  I  Two  cross  sections  through  a  highly  simplified  energy  hyper- 
surface  of  a  protein  (a)  At  constant  reaction  coordinate  (e  g  ,  B  in  panel 
b)  the  protein  possesses  two  conformations  0  and  1  (b)  At  constant 

conformational  coordinate  (e  g  ,  0  in  panel  a)  the  protein  possesses  two 
states  The  bullet  represents  the  protein-ligand  system 

at  constant  viscosity.  The  factor  expl-Pf'*gf^/PT]  can  either 
speed  up  or  slow  down  the  reaction  rate. 

As  we  will  discuss  in  more  detail  in  section  2,  a  protein  molecule 
in  a  given  state  does  not  exist  in  only  one  or  two  conformations 
but  assumes  a  large  number  of  conformational  substates  (CS)‘ 
which  most  likely  are  arranged  in  a  hierarchy Substates  arc 
defined  as  having  similar,  but  not  identical,  structure,  and  per¬ 
forming  the  same  function,  but  possibly  with  different  rates.  As 
discussed  above  for  the  simplest  case  of  two  conformations, 
pressure  can  affect  protein  dynamics  and  reactions  in  two  ways 
Within  a  given  CS,  pressure  can  change  a  reaction  rate  according 
to  eq  3  and  4;  the  change  depends  upon  the  activation  volume 
which  may  be  different  in  different  CS  F.quall>  importantly, 
however,  pressure  may  shift  the  population  from  one  substate  to 
another.  Since  the  different  CS  may  have  different  reaction  rates, 
pressure  can  change  the  overall  reaction  rate  by  a  factor  that 
cannot  be  described  by  an  activation  volume  alone. 

In  the  present  paper  we  describe  five  pressure  effects  that  are 
based  on  the  eq  1-4  and  on  the  fact  that  proteins  possess  con¬ 
formational  substates 

(i)  Well  above  a  characteristic  temperature  T,,.  defined  in 
section  2,  protein  substates  are  in  equilibrium  If  the  ratio  of 
substate  populations  can  be  monitored  as  a  function  of  pressure 
and  temperature,  then  eq  I  and  2  permit  a  determination  of  the 
differences  in  energy,  entropy,  and  volume  between  subslates  and 
thus  yield  a  detailed  thermodynamic  description  of  the  protein 
subslates  (section  4) 

(ii)  Well  below  T,,,  proteins  are  in  a  meiastable  state  The 
combined  use  of  pressure  and  temperature  proves  the  metastability 
(section  5). 

(iii)  Pressure  shifts  some  spectral  lines  These  shifts  are  caused 
by  elastic  effects  and  by  changes  in  the  populations  of  the  con¬ 
formational  subslates  (section  6) 

(iv)  Near  T^.  pressure  jump  experiments  permit  measurements 
of  the  transitions  among  various  subsiales  (section  7). 

(v)  The  binding  of  small  ligands  such  as  CO  and  Oj  to  heme 
proteins  can  be  studied  in  detail  by  using  flash  photolysis.  The 
exploration  of  this  reaction  as  a  function  of  T  and  P,  monitored 
in  the  Soret  (visible)  and  the  infrared  regions,  shows  that  pressure 
indeed  affects  reaction  rates  both  through  the  activation  volume 
and  through  shifts  in  the  substate  populations  (section  8). 

2.  Conformational  Substates  and  Motions 

2.1.  Myoglobin  Myoglobin  is  a  globular  protein,  built  of  153 
amino  acids  and  with  overall  dimensions  4.5  x  3.5  X  2.5  nm’. 
that  binds  Oj  and  CO  reversibly,'®  In  Figure  2,  we  sketch 
carbonmonoxymyoglobin  (MbCO),  with  CO  bound  to  the  heme 
iron  at  an  angle  a  with  the  heme  normal. 
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Figure  2  Schematic  view  of  myoglobin  with  CO  bound  at  the  heme  iron 
The  Fc  and  CO  are  shown  enlarged.  The  CO  dipole  makes  an  angle  a 
with  the  heme  normal 

2.2.  Stales.  Subsiales.  and  Motionsf’  Most  proteins  perform 
some  function  such  as  catalysis  or  transport  of  matter  or  charge 
They  must  consequently  possess  at  least  two  stales,  for  instance, 
Mb  and  MbOj  Proteins  in  different  states  perform  different 
functions  (Mb  binds  and  MbOj  dissociates  O;)  and  may  have 
different  structures  and  properties  As  slated  in  section  I,  each 
stale  consists  of  a  large  numlxr  of  conformational  subsiales  (CS) 
CS  have  somewhat  different  structures  but  perform  the  same 
function,  albeit  with  different  rates.  All  subsiales  of  Mb,  for 
instance,  bind  CO,  but  with  different  association  rates.  Con¬ 
formational  subslates  are  valleys  in  the  energy  hypcrsurface  of 
a  protein  in  a  particular  state. 

The  existence  of  states  and  subslates  leads  to  two  different  types 
of  motions,  equilibrium  fluctuations  (EF)  and  nonequilibnum 
transitions  called  functionally  important  notions  (FIM).’  EF 
are  motions  between  CS  in  equilibrium,  and  FIM  arc  the  none- 
quiltbrium  motions  that  lead  from  one  slate  (say  Mb)  to  another 
(MbOj)  EF  and  FIM  are  connected  by  fluctuation-dissipation 
theorems. 

2.3.  Hierarchy  of  Substates. In  MbCO.  experiments  suggest 
a  hierarchical  arrangement  of  CS  as  shown  schematically  in  Figure 

3.  where  the  conformational  energy  of  MbCO  is  plotted  as  a 
function  of  multidimensional  conformational  coordinates  (cc)  In 
principle,  listing  the  coordinates  of  all  atoms  in  a  given  protein 
molecule  describes  a  CS  completely  Figure  3  indicates  that  the 
subsiales  in  MbCO  are  arranged  in  tiers,  with  tier  0  having  the 
largest  barriers  between  energy  valleys  The  evidence  for  tiers 
in  MbCO  is  as  follows: 

Tier  0  ( eSf)  Evidence  for  the  substates  in  tier  0  comes  from 
the  stretch  frequency  vj-o  of  the  CO  bound  to  Mb.  shown  in  Figure 

4. "  The  IR  spectrum  shows  three  major  resolvable  stretch  bands 

which  wc  denote  by  Ao  (vco  =  1966  cm  ').  A,  (I'co  ~  1945  cm  '), 
and  A3  (veo  ~  1933  cm  ').  The  three  bands  correspond  to  three 
subslates;  they  have  approximately  the  same  structure  but  differ 
in  detail  as  indicated  by  the  values  of  the  till  angle  a  between  the 
bound  CO  and  the  heme  normal:  afAj)  =  15  ±  3°,  nfA,)  =  28 
±  2°,  ofAy)  =  33  ±  The  three  substates  perform  the 

same  function— binding  of  CO — but  with  different  rates:  A^,  re¬ 
binds  fastest  and  A3  slowest." 

Tier  I  ICS').  Rebinding  of  CO  to  each  CS®  is  nonexponential 
in  time,  as  proven  in  ref  1 1  and  also  shown  later  in  section  8  The 
simplest  explanation  for  the  noncxponcnlial  rebinding  is  the  ex- 
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Figure  3.  Hierarchical  arrangement  of  the  conformational  lubstates  m 
MbCO,  The  diagram  represents  a  one-dimensional  cross  section  through 
the  multidimensional  conformational  energy  surface  of  MbCO,  Three 
tiers  of  substates  (CS'.  i  =  0,  I.  2)  are  shown  as  functions  of  three 
conformational  coordinates  (cci.  /  =0.  1.2)  MbCO  can  exist  in  three 
substates  of  tier  0,  Each  of  these  can  assume  a  large  number  of  con¬ 
formations  corresponding  to  substates  of  tier  1  The  furcation  continuer, 
at  least  through  tier  2. 


Figure  4.  The  IR  spectrum  of  swMbCO  Tit  to  three  Voigtians  corre¬ 
sponding  to  three  A  bands.  The  spectrum  shown  is  for  a  sample  cooled 
to  25  K  and  then  pressurized  to  200  MPa.  pH  =  5  5  Solvent  75T 
glycerol/water  (v/v). 

istence  of  CS.'*  '*  Since  each  CS®  separately  displays  nonex¬ 
ponential  behavior,  each  CS®  must  be  furcated  into  a  large  number 
of  substates,  which  are  therefore  called  CS'.  The  nonexponential 
time  dependence  of  rebinding  proves  that  the  different  CS'  rebind 
with  different  rates. 

Tier  2  (CS^).  Originally  we  posti  lated  the  existence  of  substaics 
CS^  within  the  CS'  because  of  a  shift  of  the  near-IR  band  at  760 
nm  after  photodissociation  of  MbCO.’  Agmon'®  and  Campbell 
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et  al  '’  have  shown,  however,  that  the  major  contribution  to  the 
shift  is  due  to  a  narrowing  of  the  inhomogeneously  broadened  line 
during  the  course  of  CO  rebinding.  Nevertheless,  independent 
evidence  for  CS'  comes  from  the  Mossbauer  effect"  and  from 
pumping  experiments."  '’-^®  The  pumping  results  imply  that 
different  CS^  within  a  given  CS'  rebind  (i^O  with  different  rates. 

It  is  likely  that  substates  (CS*.  CS*)  within  CS^  exist.’  but  they 
are  not  important  for  the  present  work. 

2  4.  Slaved  Glass  Transitions  Around  300  K.  equilibrium 
lluctuations  (EFO,  EFl, ...)  occur  in  all  tiers.  As  the  temperature 
IS  lowered.  EF  in  successive  tiers  freeze  out.  As  we  have  proven 
earlier"''  and  will  show  in  sections  4  and  7,  the  exchange  between 
the  substaics  of  tier  0  stops  near  180  K  and  that  between  the 
substates  of  tier  1  near  160  K.  The  freezing  of  these  two  tiers 
is  similar  to  a  glass  transition,  but  with  one  difference:  The 
transition  temperatures  depend  on  the  solvent.  We  call  such  a 
transition  “slaved""  "  and  denote  the  transition  temperature  as 
7",,  Slaving  is  most  pronounced  for  tier  I.  where  the  transition 
in  the  protein  occurs  very  close  to  that  of  the  glass  temperature 
of  the  volvent.  Less  is  known  about  the  transitions  in  tier  2.  but 
preliminary  data  indicate  that  tier  2  freezes  out  gradually  and 
slow  exchange  may  take  place  even  at  60  K 

2.5.  Suhstates  and  Pressure  The  relevance  of  conformational 
substates  to  the  design  and  interpretation  of  pressure  experiments 
is  now  clear  for  the  binding  of  CO  or  Oj  to  Mb.  Since  different 
substates  in  each  tier  rebind  ligands  with  different  rates,  a  switch 
in  population  caused  by  pressure  affects  the  overall  binding  rate 
of  the  protein  ensemble  Pressiirc  experiments  consequently  help 
explore  protein  dynamics  in  two  different  directions.  They  provide 
information  ab*)ut  ligand  binding  and  about  the  hierarchy  of  CS. 
In  discussing  these  effects,  the  role  of  temperature  is  crucial. 
Depending  on  the  tier,  temperature,  and  reaction  or  observation 
time,  the  system  either  remains  frozen  in  a  particular  substate 
or  moves  from  CS  to  CS  In  the  frozen  ensemble,  the  reaction 
prvKceds  nonexponentially  in  time,  if  the  exchange  between 
substatcs  IS  rapid  compared  to  the  reaction,  it  shows  an  exponential 
time  dependence.'* 

3  Experimental  Methods 

Here  we  describe  the  experimental  tools  used  for  the  five  types 
of  pressure  studies  listed  m  the  ii.lioduction  The  experiments 
involve  three  types  of  measurements:  the  observation  of  the  stable 
and  mctastable  stales  of  MbCO  in  the  Soret  and  IR  regions  as 
functions  of  pressure  and  temperature,  the  observation  of  IR 
spectra  as  functions  of  time  after  a  pressure  change,  and  obser¬ 
vation  of  the  time  dependence  of  rebinding  of  CO  and  Oj  after 
flash  photolysis  as  functions  of  temperature  and  pressure  in  the 
Soret  and  IR  regions. 

l.yophilized  metmyoglobin  powders  from  sperm  whale  (sw)  and 
horse  (h)  were  purchased  from  Sigma  Chemical  Co.  Myoglobin 
solutions  were  prepared  by  dissolving  meiMb  in  75*^  (v/v)  gly¬ 
cerol/buffer  and  adding  excess  sodium  dithionite.  The  samples 
near  pH  7  were  buffered  with  O.I  M  potassium  phosphate  while 
those  near  pH  5  were  buffered  with  0.3  M  potassium  citrate.  The 
pH  of  each  IR  sample  was  measured  after  protein  was  added. 

The  pressure  cell  has  evolved  through  several  versions  over  10 
years.  Here  we  describe  the  latest  design  with  which  all  the  IR 
spectra  were  taken  (Figure  5).  The  pressure  cell  is  constructed 
of  beryllium-copper  with  two  sapphire  windows  for  Soret  and 
infrared  measurements.  The  windows  are  held  in  place  on 
mushroom  plugs  using  a  modified  "unsupported  area”  seal  which 
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Figure  5  The  high-pressure  infrared  cell  The  cell  has  been  used  at 
temperatures  between  10  and  350  K,  and  at  pressure  up  to  400  MPa  All 
parts  are  made  of  berjlhum-copper  unless  otherwise  indicated.  The  parts 
are  as  follows:  (A)  cell  body,  (B)  threaded  window  support  plug;  (C) 
stainless  steel  high-pressure  tubing  hard  soldered  to  cell  body  (0  78  mm 
outer  diameter,  0.15  mm  inner  diameter),  (D)  window  support;  (E) 
sapphire  window,  (F)  core  tube;  (G)  support  ring;  (H)  indium  packing; 
(I)  brass  antiextrusion  rings,  (J)  window  cap;  (K)  protein  sample.  (L) 
Mylar  spacer;  (M)  sapphire  windows.  (N)  sample  cell;  and  (O)  sample 
cell  bottom. 

withstands  helium  gas  pressure  up  to  400  MPa  at  temperatures 
from  40  to  350  K.  The  pressurizing  gas  enters  the  cell  through 
60  cm  of  coiled  stainless  steel  tubing  (OD  0.78  mm,  I D  0. 1 5  mm) 
enclosed  within  the  vacuum  chamber  to  minimize  heat  conduction 
along  the  tubing.  The  MbCO  sample  is  contained  between  two 
sapphire  windows  (0.25  mm  thick)  within  a  thin  doughnut-shaped 
Mylar  spacer  (0,075,  0. 1 25.  or  0. 1 75  mm).  The  sample,  sapphire 
windows,  and  Mylar  spacer  are  held  together  by  a  beryllium- 
copper  sample  holder  which  fits  inside  the  pressure  cell. 

Helium  gas  is  pressurized  by  using  a  Harwood  Engineering 
A2  5J  Intensifier  or  an  Aminco  (now  Superpressure)  compressor 
(Model  46-14060)  after  being  filtered  by  a  5-*tm  filter.  The 
pressure  is  monitored  with  a  Bourdon  gauge  accurate  to  within 
4  MPa. 

The  static  absorption  spectra  in  the  Soret  region  were  taken 
with  a  Cary- 1 4  spectrometer.  The  myoglobin  concentration  was 
about  15  fiM  and  the  path  length  was  about  7  mm. 

The  ligand  rebinding  measurements  in  the  Soret  region  were 
obtained  with  a  flash  photolysis  system.*®  The  photolysis  pulse 
was  generated  by  a  flashlamp-pumped  dye  laser  (Phase-R  DL- 
1000)  using  coumarin  6  dye  (540  nm,  100-150  mJ).  Rebinding 
was  monitored  with  light  from  a  tungsten  lamp  passed  through 
an  Interference  filter  and  collected  onto  a  photomultiplier  (RCA 
4837).  The  phot'imulliplier  output  from  7  >is  to  300  s  was  digitized 
with  a  logarithmic  time-base  digitizer.  The  pressure  cell  containing 
the  sample  was  placed  in  a  liquid  helium  storage  Dewar  (Janis), 
and  the  temperature  was  monitored  with  a  Si  diode  on  the  pressure 
cell.  The  protein  concentrations  for  the  Soret  flash  photolysis 
measurements  were  220  »iM  for  MbCO  and  300  *iM  for  MbOj; 
the  sample  path  length  was  0.5  mm. 

The  IR  spectra  were  taken  on  a  Mattson  Sirius  100  FTIR 
spectrometer  at  2  cm*'  resolution.  Low  temperatures  were 
maintained  with  a  closed-cycle  helium  refrigerator  (CTl  Model 
22C)  and  a  Lake  Shore  Cryogenics  temperature  controller  (Model 
93C).  A  silicon  diode  mounted  on  the  sample  cell  monitored  the 
temperature.  Photoly.ls  for  the  IR  rebinding  experiments  was 
achieved  with  an  argon  ion  laser  operating  at  514  nm  with  an 
intensity  of  0.4  W/cm^  and  a  beam  diameter  of  about  3  mm.  The 
pH  6.6  samples  were  15  mM  MbCO,  and  were  photolyzed  for 
20  s;  the  pH  5.5  samples  were  5  mM  MbCO  and  were  photolyzed 
for  10  s.  All  samples  used  for  flash  photolysis  were  0,075  mm 
thick.  Absolute  absorbance  spectra  (as  in  Figure  9)  were  obtained 
by  subtracting  either  a  solvent  background  or  a  solvent  deoxy- 


T(K) 


Figure  6  Logarithm  of  the  ratio  of  areas,  Ao/A|,  of  CO  stretch  bands, 
at  pH  5  5  (a)  and  pH  6  6  (b).  as  a  function  of  lOOO/T  for  different 
pressures  Sample.  swMbCO  in  75%  glyuciol/ water  (v/v) 

myoglobin  background  from  the  MbCO  spectrum  Sperm  whale 
myoglobin  was  used  for  all  IR  samples. 

The  CO  stretch  bands  were  fit  in  a  two-step  process.  First,  the 
base  line  was  fit  to  a  cubic  polynomial  by  a  standard  linear 
least-squares  method.  Then,  the  bands  were  fit  to  Voigiian  line 
shapes  (Gaussian  superpositions  of  Lorentzians)  by  using  a  finite 
difference  Levenberg-Marquardt  nonlinear  least-squares  algor¬ 
ithm.  The  spectra  were  fit  to  four  Voiglians  to  account  for  the 
asymmetry  of  Ao  at  pH  5.5  and  of  A,  at  pH  6.6. 

4.  Equilibrium  Properties 

The  behavior  of  the  substates  of  tier  0  as  a  function  of  Rand 
T  is  shown  in  Figure  6.  w  here  the  ratios  of  the  areas  of  the  CO 
stretch  bands  of  swMbCO.  for  two  pH  values,  are  plotted  vs 
1000/r  for  different  pressures.  Above  about  330  K  MbCO 
unfolds,  and  the  entire  CO  stretch  spectrum  changes.  Between 
about  330  and  200  K,  the  three  A  subslates  are  in  equilibrium, 
in  the  transition  region  between  about  2(X)  and  180  K,  exchange 
between  the  A  subslates  becomes  measurably  slow.  Below  about 
180  K,  the  excha.nge  is  immeasurably  slow  and  each  protein 
molecule  remains  frozen  into  a  particular  CS®;  the  protein  has 
become  glasslike.  We  use  the  terms  glassy  or  glasslike  to  describe 
a  frozen  and  metastable  complex  system.  In  the  present  section, 
we  treat  the  equilibrium  region. 

The  equilibrium  region  of  Ao/A,  shows  two  van't  Hoff-like 
parts,  with  a  minimum  near  300  K  We  have  shown  that  for 
constant  pressure  this  behavior  follows  from  the  thermodynamic 
properties  of  proteins. Here  we  generalize  the  treatment  to 
include  pressure  and  shuw  that  the  data  in  Figure  6  lead  to  a 
determination  of  the  relative  thermodynamic  parameters  of  the 
A  substates. 

We  assume  that  the  A  substates  have  equal  extinction  coef¬ 
ficients.^’  The  ratio  of  the  populations  of  the  substates  in  eq  1 
is  then  equal  to  the  ratio  of  the  band  areas  If  the  thermodynamic 
parameters  E,  V.  and  5  are  taken  to  be  independent  of  T.  the  data 
in  Figure  6  cannot  be  fit  with  only  the  two  substates  Ao  and  A,. 
Experiments  show,  however,  that  the  specific  heat  of  globular 
proteins  is  essentially  linear  from  200  to  320  K.’*  The  entropy 
can  consequently  be  approximated  hy  S(T)  =  S(0)  +  sT.  We 
generalize  this  expression  to  also  include  a  term  linear  in  pressure 
S{T.P)  =  S(0)  +  sT  +  iP  (5) 
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Table  I:  Paramcten  for  Ibc  Population  Ration 

A«/A,  at  pH  6.6  and  5.5* 


pH 


parameter 

6  6 

5.5 

A£(0).  kj  mol-' 

-68 

-22 

AKO).  cm’  mol"' 

-45 

-42 

AA(0),  J  mol-'  K  ' 

-510 

-150 

Ac,  cm’  mo|-'  K-' 

-Oil 

-0.13 

Ar,  J  mol-'  K-’ 

+  17 

+0.44 

•Data  and  fit  sho'vn  in  Figure  6.  Errors  are  ±25%.  Sample: 
swMbCO  in  75%  glycerol/water  (v/v). 

This  expansion  is  justified  by  noting  that  the  Maxwell  relation 
(dS/dP)T=  -(dl'/dr),  yields 

y{T.P)  =  y(0)  -  vT  (6) 

With 

f  =  ~3y(0)  (7) 


eq  6  is  the  standard  relation  F'(r)  =  F'(O)  (1  +  371.  where  3  is 
the  coefficient  of  thermal  expansion  and  where  the  compressibility 
has  been  neglected.  From  eq  5,  standard  thermodynamics  yields 

E(T.P)  =  £(0)  +  lPT  +  %s7^  (8) 

The  difference  in  Gibbs  free  energy  between  substates  Ao  and  A, 
becomes 

NG(T.P)  3  Gq-C,  = 

A£(0)  +  P  AF'(O)  -  r  A5(0)  -  \i  PT  ~  '/^\sT^  (9) 

where,  for  instance,  A£(0)  =  £o(0)  -  £((0),  and  £o(0)  is  the 
internal  energy  of  substale  Ag  at  T  =  0,  £  =  0. 

Fits  to  eq  I  and  9  yield  the  solid  lines  in  Figure  6  and  the 
parameters  in  Table  I.  These  values,  together  with  eq  5,  6,  and 
8,  give  the  differences  in  entropy  (in  units  of  /?),  energy,  and 
volume  between  substates  as  a  function  of  temperature  shown  in 
Figure  7. 

In  order  to  discuss  the  results  in  Table  I  and  Figure  7  we  first 
note  that  the  volume  differences  are  given  in  cm’/mol,  with  1 
cmVnaol  =  1.67  X  10"’  nmVmolecule  =  1.67  A’/ molecule  The 
volume  F'y,,  of  a  Mb  molecule  is  about  40000  A^  The  pocket 
volume  is  about  200  A’.  The  ratios  IV are  of  the  order  of 
10"’,  indicating  that  the  overall  volumes  of  the  three  substates 
are  nearly  the  same.  The  volume  of  Ag  is  smallest,  but  the 
difference  A  F  decreases  with  increasing  temperature.  Below  about 
300  K,  Ag  also  has  the  smallest  entropy  and  the  lowest  internal 
energy.  This  observation  suggests  that  Ag  may  be  stabilized  by 
additional  hydrogen  bonds  at  low  T  and  that  the  number  of  these 
bonds  decreases  with  increasing  T.  X-ray  diffraction  studies  at 
different  pFI  and  pressures”-’*  and  calorimetric  studies  at  different 
pFI  are  needed  to  explore  the  substates  in  more  detail. 

A  last  remark  here  concerns  the  value  of  v.  With  eq  7  we  write 
Ap  =  -A(/3K(0)1  -/3AF(0)  -  F(0)A/j  (10) 

Using  AF'(0)  from  Table  I  and  taking  the  volume  thermal  ex¬ 
pansion  coefficient  3  fiom  X-ray  data”  as  3.5  X  10"*  K"'  gives 
for  the  first  term  in  eq  10  the  value  -FO.Ol  cmV(mol  K).  Since 
Ap  for  Ag  is  negative,  the  main  contribution  to  Ap  must  come  from 
the  second  term,  suggesting  that  A/J/fl  =  Ap/^FfO)  =  Ap//3F^,^ 
=  0.005. 

S.  The  Mctastablc  Glass  Slate 

Below  about  160  K,  the  exchange  among  substates  of  tier  0 
and  tier  I  is  frozen,  and  the  protein  is  glasslike. Two 
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Figure  7.  Differences  m  entropy  (in  units  of  R).  energy  (kJ/mol),  and 
volume  (cm’/mol)  between  the  substates  Ag  and  A,  at  pH  6.6  and  pH 
5.5  at  zero  pressure  as  a  function  of  temperature.  The  differences  are 
derived  from  the  fits  in  Figure  6. 
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Figure  8.  The  protein  sample  is  taken  from  an  initial  slate  {T..P.)  to  a 
final  state  {Tj.P,)  on  different  pathways,  for  instance  0  —  F  —  FP  (solid 
path)  and  0  -»  P  —  PF  (dashed  path).  In  equilibrium,  the  protein 
properties  in  the  two  slates  must  be  identical.  Metastability  is  proven 
if  different  properties  are  found  al  FP  and  PF. 

outstanding  characteristics  of  the  glassy  (amorphous)  state  are 
a  rugged  energy  landscape  and  metastability.  The  rugged  energy 
landscape,  valleys  within  valleys  within  valleys,  is  shown  in  Figure 
3.  Metastability  implies  that  the  state  of  the  system  below  the 
glass  temperature  depends  on  its  history  and  that  the  system  is 
trapped  in  a  substate  with  an  energy  above  the  lowest  energy. 
History  dependence  can  be  shown  by  moving  in  the  TP  plane 
(Figure  8)  from  a  point  above  the  glass  transition  to  a  point  below 
by  two  different  pathways.’*  In  equilibrium,  the  properties 
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Figure  9.  The  MbCO  stretch  spectrum  determined  at  the  points  F,  FP, 
PF.  PFR  in  the  TP  plane  of  Figure  8,  with  T,  =  225  K,  T,  =  100  K,  P, 
=  0.1  MPa,  and  Pf  =  200  MPa  In  equilibrium,  the  pairs  F/PFR,  and 
FP/PF  should  yield  identical  spectra  The  different  spectra  prove  me- 
tastabilily.  Sample  swMbCO  in  75%  glycerol/watcr  (v/v)  at  pH  5  5. 

of  the  system  are  independent  of  the  path  (history)  For  most 
of  our  experiments  we  select  the  two  pathways  shown  in  Figure 
8.  On  the  solid  pathway,  we  freeze  the  sample  (F)  and  then 
pressurize  (FP).  On  the  dashed  pathway,  we  first  pressurize  (P). 
then  freeze  (PF),  and  finally  release  the  pressure  (PFR).  In 
equilibrium.  F  and  PFR  or  FP  and  PF  should  yield  identical 
protein  properties.  Figure  9  shows  the  results  of  such  an  ex¬ 
periment,  where  the  CO  stretch  spectrum  of  swMbCO  is  measured 
at  the  four  points  F,  FP,  PF,  PFR  in  the  TP  plane,  with  T,  =  225 
K.  T,  =  100  K,  P,  =  0.1  MPa,  P,  =  200  MPa.  The  infrared 
spectra  obtained  via  the  two  pathways  differ  considerably. 
Metastabilily  is  proven. 

6.  Pressure  Effects  on  Spectral  Lines 

6./.  Elastic  and  Conformational  Pressure  Effects.  In  crystals, 
the  dominant  effect  of  a  moderate  pressure  on  the  optical  spectra 
is  a  peak  shift,  usually  to  lower  wavenumber  (red  shift).  The 
pressure  compresses  the  crystal  and  decreases  the  distances  be¬ 
tween  the  atoms,  thereby  changing  the  total  energy.  The  ar¬ 
rangement  (topology)  of  the  atoms  remains  unchanged.  Since 
this  effect  involves  only  an  elastic  deformation  of  the  crystal,  it 
occurs  rapidly  at  any  temperature.  At  the  pressures  of  interest 
here  {P  <  200  MPa),  the  wave  number  shift  is  linear  in  P 

S‘'.  =  c,P  (ID 

and  the  coefficient  c,  is  proportional  to  the  volume  difference 
between  the  ground  and  the  excited  state  of  the  optical  transition.’* 

in  proteins,  the  situation  is  more  complex  and  two  different 
types  of  pressure  effects  occur,  elastic  and  conformational 
(plastic).’*  The  elastic  effect  is  similar  to  the  one  in  crystals; 
The  protein  remains  in  the  same  substate  and  does  not  change 
conformation,  and  the  arrangement  (topology)  of  the  atoms  re¬ 
mains  unchanged.  The  spectral  shift  is  due  to  an  increase  in  the 
density  and  hence  energy  of  the  protein.  The  conformational  effect 
is  characteristic  of  proteins,  but  is  similar  to  the  effect  of  pressure 
on  thermal  transitions  discussed  by  Drickamer  and  co-workers.” 
In  conformational  effects,  pressure  changes  the  local  or  global 
arrangement  of  the  atoms;  the  protein  moves  from  one  confor¬ 
mational  substate  to  another.  Fjtch  conformational  substate  may 
give  rise  to  a  line  with  a  slightly  different  peak  wavenumber,  and 
the  overall  band  may  be  a  Gaussian  superposition  of  Lorentzians 
(Voigtian).  Thus,  spectral  lines  in  proteins  are  inhomogeneously 
broadened;'*  ”-’*”  a  pressure-induced  shift  in  the  substate 
populations  not  only  shifts  the  line  but  also  changes  its  shape. 


(35)  Slichler,  C.  P.;  Drickamer,  H,  G.  Phys.  Rev  B  19*0,  22,  4097. 

(36)  Noguti,  T.;  Go,  N.  Proteins  1989,  5,  97. 

(37)  Drickamer.  H.  G.;  Frank.  C.  W.;  Slichter.  C.  P  Proc  Nall.  Acad. 
Sci.  USA  1972,  69,  933. 

(38)  Agmon,  N.;  HopField,  J  J  J  Chem  Phys  19*3,  79.  2042 

(39)  Ormos,  P.;  Anuri,  A.;  Braunstein,  D.;  Cowen,  B.  R.;  Frauenfelder. 
H.;  Hong,  M.  K.;  Iben,  I.  E.  T.;  Sauke,  T.  B.;  Steinbach,  P.  J.;  Young.  R  D 
Biophys.  J.,  in  press. 


Solvent 

Pretluft 

Proie*n 

P 


LiQOrtd 

Pressure  P|,|, 
Botind  Wottr 


Figure  10  Pressure  screening  An  externally  applied  pressure  may 
not  be  fully  transmitted  to  the  protein.  The  pressure  seen  by  the  protein 
is  given  by  P  =  qPo,  The  ligand  bound  at  the  heme  may  lx  additionally 
shielded  by  the  heme  pocket  so  P^^.  the  pressure  the  ligand  feels,  may 
differ  from  P. 


Elastic  and  conformational  band  shifts  can  be  distinguished  with 
the  technique  depicted  in  Figure  8;  At  the  temperature  T,  (>T„) 
both  the  elastic  and  the  conformational  shifts  occur.  At  T,  («tJ), 
only  the  elastic  shift  takes  place.  To  discuss  these  shifts,  we 
consider  the  simplest  case  and  assume  two  conformational  sub¬ 
states  (Figure  la),  with  peak  wavenumbers  vo  and  v,  and  elastic 
shifts  6ff)  and  Si',.  We  further  assume  the  normalized  populations 
woiT)  and  ycfP)  of  the  two  substates  to  be  functions  of  pressure. 
The  average  wavenumbers  in  the  four  pressure  states  of  Figure 
8  then  become 

I'F  =  Wo(0)vo  +  W|(0)l-| 

fFP  =  ''b(0)(>'o  +  ^I'o)  *'i(0)(r'i  +  ^<'i) 

Ppr  =  *v£|(R)(vo  +  ivo)  +  “  i(f*)(Pi  +  i<'i) 

*'PFR  =  '*o(^)‘'o  +  ‘*i(^)>'i  (12) 

From  these  four  observables  we  form,  with  wo(/’)  +  h'i(/’)  =  1, 
three  independent  differences 

A|  =  Vff.  -  »'f  =  »'o(0)^‘'o  +  w'|(0)6p| 

A;  =  I'PF  -  VpFR  =  W-ofRISl'o  +  H|(/’)6V| 

•^3  ^  '*'PFR  "  *'F  ~  |w'|(f’)  ~  '*|(0)]|V|  -  V[)]  (13) 

The  first  two  differences  are  elastic,  and  they  vanish  if  no  elastic 
shift  occurs;  i.e.,  if  ipp  =  6pi  =  0.  The  third  difference  is  con¬ 
formational  and  is  nonzero  only  if  the  pressure  changes  the 
population  of  the  substates  and  the  wavenumbers  of  the  substates 
differ.  Finally,  a  nonvanishing  difference  between  A;  and  A2 

•^2-Ai  =  |h-i(/’)- w,(0)][ap, -iPp)  (14) 

proves  that  even  the  elastic  shifts  in  the  different  substates  differ. 

6.2.  Pressure  Screening.  In  a  crystal,  the  pressure  inside  is 
the  same  as  the  externally  applied  hydrostatic  pressure  P„,.  In 
a  protein,  the  situation  is  more  intricate  as  indicated  in  Figure 
10.  Two  effects  must  be  taken  into  account;  In  the  glassy  state 
below  the  glass  temperature,  T,,  of  the  solvent,  the  solvent  may 
not  transmit  the  entire  pressure  to  the  protein.  At  all  temperatures, 
the  ligand  within  the  pocket  may  be  shielded  even  from  the  protein 
pressure  P  by  the  pocket. 

In  a  simple  static  model,  where  the  coefficients  of  thermal 
expansion  are  neglected,  the  effective  pressure  at  the  protein  can 
he  calculated.*®  The  resulting  effective  pressures  in  the  four  states 
below  T,g  (Figure  8)  are 

Pf~0.  Pfy=^qP^^,,  RpF  "  ^Ml-  ^PFR  =  (1  “  9)^eil  (15) 
where 

1  +  (4/3)u,<c, 

<7  =  7-7777^7 -  (16) 

1  +  (4/3)m,«p 

Here  p,  is  the  shear  modulus  of  the  solvent,  and  x,  and  Xp  are  the 
compressibilities  of  the  solvent  and  the  protein.  Protein  com- 


(40)  Sorensen,  L  B.  Pti.D.  Tliesis,  University  of  Illinois  at  Urbina- 
Champaign,  1980 
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Figure  II  The  Sorei  band  for  different  pressure  states.  Sample; 
swMbCO  in  75%  glycerol/water  solvent  by  volume,  pH  7  T  =  140  K, 
P  =  190  MPa. 

TABLE  II:  Peak  Position  of  the  Soret  Band  for  swMb  in  the 
Pressure  Slates  Shown  in  Figure  8* 

peak  position,  nm 


protein 

T.  K 

solvent 

F 

FP 

PF 

PFR 

MbCO 

140 

water 

423.7 

424.4 

424.4 

423  5 

MbCO 

140 

75%  G 

423.4 

423  8 

424.5 

423.8 

MbCO 

140 

97%  G 

423  5 

424.1 

424.4 

423.5 

Met  Mb 

120 

75%  G 

4100 

410.7 

411  3 

410.2 

"Pressure  at  P  is  190  MPa.  G  denotes  glycerol.  %  by  volume.  50 
mM  phosphate  buffer,  pH  7.  Points  F,  FP,  PF,  and  PFR  are  explained 
in  Figure  8. 

pressibilites  have  been  measured  by  using  sound  velociiy*'-^^  and 
fluorescence^^  techniques.  The  sound  velocity  technique  measures 
an  average  value,  the  fluorescence  experiment,  a  change  in  distance 
between  the  heme  and  a  tryptophan  residue.  Since  the  distance 
can  also  be  changed  by  a  pressure-induced  change  in  substates, 
and  since  in  eq  16  refers  to  an  average  protein  value,  we  use 
the  results  from  the  sound  velocity  experiments.  Some  of  the  data 
in  the  next  subsection  refer  to  ice  as  the  solvent.  The  relevant 
coefficients  for  ice  are*^  m,  =  3.8  GPa  and  x,  =  O.1 1  GPa"';  the 
compressibility  of  Mb  lies  between  0.06  and  0.1 2  GPa"'.  For  ice 
as  the  solvent,  q  lies  between  1.2  and  0.97.  For  the  following 
discussions  we  assume  q  =  \. 

We  will  see  below  that  there  exists  evidence  for  additional 
screening  of  the  CO  molecule  bound  to  the  heme  iron.  The 
screening  depends  markedly  on  substate.  To  interpret  the 
screening.  X-ray  structures  for  the  different  CS®  will  be  required. 

Differences  in  the  thermal  expansion  of  the  protein  and  the 
solvent  can  also  induce  a  pressure  on  the  protein.”  If  the  thermal 
expansion  coefficient  of  the  protein  is  smaller  than  that  of  the 
solvent,  a  temperature  change  may  lead  to  tension,  cancelling  part 
of  the  external  pressure.  The  experiments  that  will  be  discussed 
in  section  8  indicate  that  this  effect  is  not  pronounced. 

6.3.  Pressure  Shift  of  the  Soret  Band.  A  pressure-induced 
red  shift  of  the  Soret  band  in  heme  proteins  has  been  observed 
by  a  number  of  groups.^^’  Since  these  experiments  were  per¬ 
formed  in  a  small  temperature  range,  they  do  not  provide  in¬ 
formation  on  the  separation  of  the  shift  into  elastic  and  confor¬ 
mational  components. 

We  have  measured  the  Soret  band  of  swMb  in  a  number  of 
solvents  at  the  four  points  F,  FP,  PF,  and  PFR  of  Figure  8.  One 
set  of  data  is  shown  in  Figure  1 1,  and  the  peak  wavelengths  in 
a  few  solvents  are  listed  in  Table  II.  For  example,  the  values 
for  swMbCO  in  a  75%  glycerol-water  solvent  (Figure  1 1 )  give 


(41)  G«kko.  K.;  Noguchi,  H.  J  Phys.  Chem.  197».  S3.  2706. 

(42)  Gavish,  8.;  Gration,  E.;  Hardy,  C.  J.  Proc.  Nall.  .lead.  Sci.  USA 
IW3,  SO.  750. 

(4))  Marden,  M.  C.:  Hui  Bon  Hot.  G.:  Stetzkowiki-Marden,  F.  Biophys. 
J  I9M.  49,  619. 

(44)  Hobbt,  P.  V,  let  Phystcr  Oxford:  London,  1974. 

(45)  Ogunmoli,  G.  8.;  Zipp,  A.;  Chen,  F.;  Kauzmann,  W.  Proc.  Nall. 
Acad.  Sel  USA  1*77,  74.  1. 

(46)  Gibion,  Q  H.:  Carey,  F.  G.  J  Biol.  Chem.  1977,  2S2,  4098. 
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TABLE  IH;  Wavenumber  Differences  for  the  Stretch  Bands  A*  and 
A,  of  swMbCO" 


elastic  conformational 
shifts,  cm  '  shifts,  cm"' 


substate 

pH 

T.  K 

A, 

Ay 

A, 

Ay  -  A, 

Ao 

5.5 

25-100* 

-16 

-0.7 

-0.2 

+09 

6.6 

50 

-1.3 

-0  5 

-1.4 

+0.8 

A, 

5.5 

25-100* 

0 

66 

50 

0 

-0.1 

+0.3 

-0.1 

‘  Differences 

are 

given  in  cm 

'.  errors  ±0  2 

cm'*. 

Pressure 

MPa.  75%  glycerol/walcr  (v/v).  Potassium  phosphate  buffer.  The 
points  F,  FP,  PF,  and  PFR  are  explained  in  Figure  8  ‘Average  values. 

the  following  with  eq  1  3  and  14  and  with  -A\/X  =  Se/v:  S,  = 
-20  cm"';  Ay  =  -40  cm"'.  A)  =  -20  cm"'.  Ay  -  A,  =  -20  cm"'. 
The  nonzero  values  of  A,  and  Ay  prove  that  elastic  shifts  occur 
in  the  conformations  between  P  =  0  and  P  =  190  MPa.  The  third 
difference  implies  that  pressure  changes  the  populations  and  that 
the  peak  wavenumbers  of  the  substates  are  different.  The  finite 
value  of  Ay  -  A,  demonstrates  that  the  elastic  shifts  in  various 
substates  are  different 

The  shifts  for  .VfbCO  in  75%  glycerol/water  have  been  mea¬ 
sured  in  the  temperature  range  from  60  to  160  K  and  were  found 
to  be  independent  of  T.  This  independence  implies  that  the 
pressure  effects  caused  by  thermal  expansion  are  small. 

The  pressure  effect  on  the  Soret  band  does  not  determine  in 
which  tier  of  substates  the  pressure  effect  occurs.  This  information 
can  be  obtained  from  pressure  studies  of  the  CO  stretch  bands 
in  MbCO 

6.4.  Pressure  Effects  on  the  CO  Stretch  Bands.  Studies  of 
the  stretch  vibrations  in  condensed  systems  by  Drickamer  and 
collaborators'** show  that  hqq  red  shifts  with  increasing  pressure; 
the  coefficient  Cj  in  eq  1 1  is  of  the  order  of  3  X  10"^  cm"'  MPa"'. 
We  have  studied  the  effect  of  pressure  on  the  stretch  bands  Ag 
and  A|  of  swMbCO.  .An  example  is  shown  in  Figure  9.  Apart 
from  the  change  in  population  already  discussed  in  section  4, 
pressure  shifts  the  stretch  bands  and  changes  their  shape.  The 
shifts  and  shape  changes  depend  on  the  path  in  the  TP  plane.  In 
Table  III  we  list  the  wavenumber  differences  as  defined  in  eq  13 
and  14  for  the  two  bands  Ag  and  A|  at  two  different  pH  values. 

The  two  substates  Ag  and  A|  are  surprisingly  different.  A, 
shows  very  small  or  vanishing  elastic  and  conformational  shifts. 
Ag.  in  contrast,  displays  both  elastic  and  conformational  shifts. 
Even  for  Ag,  however,  the  coefficient  Cj  in  eq  1 1,  c,  =  5  X  10"’ 
cm"'  MPa"’,  is  considerably  smaller  than  for  carbonyls  in  liquids. 
The  smallness  of  the  elastic  shifts  and  the  different  behavior  of 
Ag  and  A,  may  be  caused  by  different  pressure  screening  or  by 
a  different  dielectric  environment  within  the  heme  pocket.  The 
Kirkwood-Bauer-Magat  relation’®  gives  for  the  elastic  shift 

f  -  1 

—  =  const  X  — —  (17) 

v  2t  -F  1 

where  « is  the  dielectric  coefficient  of  the  medium  surrounding 
the  CO  dipole.  Different  c  may  account  for  the  magnitudes  and 
differences  of  the  elastic  shifts. 

The  conformational  pressure  shift  of  Ag  implies  that  the  sub¬ 
states  of  tier  1  within  a  given  CS°  have  properties  different  enough 
to  produce  a  conformational  shift.  This  conclusion  is  supported 
by  three  more  observations:  The  substate  Ag  broadens  under 
pressure  as  shown  in  Figure  9,  the  tilt  angle  a  between  the  CO 
axis  and  the  heme  normal  In  a  given  CS®  is  distributed  about  its 
average  value,'’  and  different  parts  within  a  given  CO  stretch  band 
rebind  CO  after  photodissociation  with  different  rates.’* 

7.  Pressure  Studies  of  Protein  Dynamics 

So  far  we  have  discussed  the  application  of  pressure  to  the  study 
of  protein  structure,  equilibrium  behavior,  and  metastability. 
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Figure  12.  The  behavior  of  an  observable  0{P:l)  after  a  rapid  pressure 
change  at  time  i  =  0.  An  essentially  instantaneous  elastic  relaxation  is 
followed  by  a  slower  conformational  relaxation. 

Pressure  is  also  a  good  tool  wi  It  which  to  explore  the  dynamics 
of  proteins 

7.1.  The  Pressure  Jump  Technique.^'  In  the  transition  region 
shown  in  Figure  8  between  about  150  and  210  K,  the  relaxation 
phenomena  are  slow  enough  to  permit  the  determination  of  the 
relaxation  functions  with  a  Fourier-transform  infrared  spectrom¬ 
eter.*' 

The  idea  of  a  pressure  jump  experiment  is  simple.  Consider 
an  observable  0(P,T,l)  that  depends  on  pressure,  with  an  equi¬ 
librium  value  0{P.T.<^)  =  0(P,T).  The  system  is  first  equilibrated 
at  a  given  temperature  T  with  pressure  P,.  At  time  /  =  0,  the 
pressure  is  changed  rapidly  to  a  new  value  Pf  and  the  change  of 
the  observable  OiP.T.i)  as  a  function  of  time  after  the  pressure 
change  is  measured.  The  general  behavior  of  0{P.T\i)  after  such 
a  pressure  jump  is  shown  schematically  in  Figure  1 2.  Initially 
0{P,T\i)  changes  very  rapidly  from  0{P,,T)  to  a  value  OiPf.T.O*). 
We  identify  this  change  as  elastic  relaxation.^  In  a  second  phase, 
OiP.T.i)  relaxes  with  measurable  rate  toward  the  new  equilibrium 
value  OiPf.T'^).  We  interpret  the  second  change  as  the  con- 
formational  relaxation  and  characterize  it  at  the  temperature  T 
by  the  relaxation  function  4>(T;r), 


*iTt)  = 


OiP,.T.t)  -  OiP^.T.’^) 
OiPf.TO*)  -  OiPf.T'^) 


The  experimental  determination  of  4>{T;t)  is  made  difficult  by 
the  nonexponential  character  of  the  conformational  relaxation. 
After  the  pressure  jump,  the  experiment  must  be  continued  long 
enough  so  that  the  final  value  OiPf.T<^)  is  sufficiently  well-known. 

7.2.  Time  and  Temperature  Dependence  of  <l>.  The  relaxation 
functions  observed  in  complex  systems  are  usually  not  exponential 
in  time  and  cannot  be  characterized  by  a  single  rate  coefficient.”'^ 
It  is  customary  to  approximate  their  time  dependence  either  by 
a  stretched  exponential 

4>(r;r)  =  exp[-(Mr)fn  (19) 

or  by  a  power  law'^ 

*iTt)  =  [\  -i- kATit]-"  (20) 

The  parameters  k,,  0,  and  n  in  eq  19  and  20  are  usually  functions 
of  T.  The  rate  parameter  k^iT)  follows  an  Arrhenius  relation  only 
over  small  intervals  in  T.  In  glasses,  k,(T)  can  be  described  over 
more  than  10  orders  of  magnitude  either  by  the  Vogel-Tam- 
mann-Fulcher  equation” 

k,iT)  =  Aytxp[-Ey/R{T-Ty)]  (21) 

or  by  the  relation”” 

k,{T)  =  A^expl-iEt/RD^]  (22) 
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Figure  13  The  peak  wavenumber  i<o(r)  of  the  Ao  band  plotted  versus  log 
/  in  a  pressure-release  experiment  Top  Pressure  plotted  versus  log  r. 
Bottom:  (b)  150-170  K.  The  elastic  shift  (0,5  cm"')  is  shown  together 
with  the  prerelease  wavenumber  (a)  175-205  K  Sample.  swMbCO 
in  75“*  glycerol/water  at  pH  6.6 

Equation  22  contains  only  two  free  parameters  and  consequently 
yields  better  determined  parameters  than  eq  21  if  data  are 
available  only  over  small  temperature  intervals. 

7.3.  Relaxation  Processes  in  suMbCO.  The  results  presented 
in  Figures  6  and  9  provide  us  with  three  observables  O(P.t)  for 
studying  the  relaxation  processes  in  MbCO.  the  relative  areas, 
the  peak  wavenumbers,  and  the  widths  of  the  CO  stretch  bands 
A,.  The  information  in  seciion  2  permits  a  logical  connection 
between  these  observables  and  the  tiers  of  substates.  The  three 
CO  stretch  bands  label  the  three  substates  of  tier  0.  Exchange 
between  these  CS”  after  a  perturbation,  for  instance  a  pressure 
jump,  hence  characterizes  FIM  0,  the  nonequilibrium  motion 
among  CS°.  The  peak  wavenumber  and  width  of  an  individual 
CO  stretch  band  characterizes  the  subslates  of  tier  1  within  a  given 
CS°.  We  therefore  assume  that,  in  the  absence  of  exchange  with 
another  CS",  the  motions  of  these  observables  after  a  perturbation 
label  FIM  1,  the  nonequilibrium  motion  among  the  CS'  of  a  given 
CS®. 

We  have  already  published  a  set  of  relaxation  data  obtained 
with  the  pressure-release  technique.^'  These  experiments  were 
performed  with  swMbCO  in  75%  glycerol/water  (v/v)  at  pH  = 
7  by  applying  ICX)  MPa  at  240  K,  cooling  the  sample  to  the  desired 
tem;>erature  in  the  transition  region  between  160  and  205  K, 
measuring  the  IR  spectrum  for  about  300  s,  and  releasing  the 
pressure  to  7  MPa  in  a  few  seconds.  IR  spectra  were  then  taken 
at  approximately  exponentially  increasing  times  from  10  to  more 
than  10‘  s  The  relaxation  functions  ♦off)  and  ♦)(()  were  de¬ 
termined  by  eq  18  where  O  is  the  area  of  band  Ao  for  ♦o,  and 
the  peak  wavenumber  or  width  of  Ao  for  ♦,.  The  resulting  re¬ 
laxation  functions  are  shown  in  ref  21  and  reveal  two  charac¬ 
teristics  that  are  also  seen  in  relaxation  processes  in  glasses:  The 
relaxation  functions  are  nonexponential  in  time,  and  the  relaxation 
rates  differ  by  orders  of  magnitude  in  a  narrow  temperature 
interval  suggesting  a  non-Arrhenius  temperature  behavior.  The 
time  dependence  of  ♦,(/)  was  fit  by  using  a  power  law,  eq  20,  which 
describes  the  data  better  than  a  stretched  exponential,  eq  19. 

A  new  set  of  pressure  release  data  on  MbCO  in  75%  glycer¬ 
ol/water  at  pH  5  5  and  6,6  again  shows  the  two  relaxation  pro¬ 
cesses  FIM  0  and  FIM  1  but  surprisingly  also  gives  evidence  for 
a  third  relaxation  process  near  180  K.  The  data  were  taken  by 
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Figure  14.  The  areas  of  the  bands  Ao  and  A,  as  a  function  of  time  after 
a  pressure  release  from  200  to  25  MPa.  Temperature,  180-205  K. 
Sample;  swMbCO  in  glycerol/ water  at  pH  6  6 


applying  200  MPa  at  225  K,  cooling  the  sample  to  the  desired 
temperature,  measuring  the  IR  spectrum  for  about  300  s,  and 
releasing  the  pressure  to  25  MPa  in  a  few  seconds.  Infrared 
spectra  were  again  taken  at  approximately  exponentially  increasing 
times  from  about  10  to  10*  s.  Figure  13  shows  the  peak  wave- 
number  Vo  for  band  Ao  at  pH  6.6  from  1 50  to  205  K  plotted  versus 
log  I.  Figure  14  shows  the  areas  of  Ao  and  A|  as  a  function  of 
time  after  pressure  release.  The  data  in  these  figures  suggest  the 
following  sequence  of  relaxations  in  the  transition  region  near  Tj,: 

„  ,  -  «>«<'«  ,  .  FIM  I  „  _  Fist  X  ^  FIM  0 

1965.7  - -  1966.2  - -  1967.0  - -  1966.2  - - 

1966.2  cm-' 


After  the  pressure  release,  a  rapid  elastic  shift  changes  the  peak 
position  of  Ao  from  1965.7  to  about  1966.2  cm"'.  Below  about 
150  K,  no  further  shift  occurs.  Above  about  150  K,  the  peak 
position  of  Ao  does  not  stop  at  1966.2  cm"'  but  continues  to  shift 
to  the  blue  with  a  nonexponential  time  dependence.  Above  180 
K,  FIM  1  is  so  fast  that  relaxation  appears  to  start  at  1967.0  cm"’. 
Surprisingly,  however  1967.0  cm"'  does  not  correspond  to  the 
relaxed  pressure  state  F,  1966.2  cm"'.  Indeed,  a  new  relaxation 
process  occurs  which  leads  from  1967.0  to  1966.2  cm"'.  Since 
we  do  not  yet  have  a  satisfactory  interpretation  of  this  relaxation 
process,  we  call  it  FIM  X.  Above  180  K,  FIM  0  sets  in;  it  is  the 
transition  from  substate  Ao  to  A|  measured  by  the  change  in  area 
of  A.  Above  about  210  K,  FIM  0  becomes  too  fast  to  be  observed 
with  our  system. 

In  addition  to  the  three  relaxation  processes  FIM  0,  1,  and  X 
at  least  one  other  similar  process  occurs  near  180  K:  the  exchange 
between  the  substates  A|  and  Ay.^'  We  have  not  yet  investigated 
this  exchange  in  detail;  its  study  will  require  pH  greater  than  7. 

To  evaluate  the  relaxation  data,  we  first  calculate  the  relaxation 
functions  ♦ft),  eq  18.  The  observable  O  in  eq  18  is  the  peak 
frequency  or  width  of  Aq  for  FIM  1,  the  peak  frequency  of  Ao 
for  FIM  X,  and  the  area  of  A©  for  FIM  0.  The  three  relaxation 
functions,  displayed  in  Figure  IS,  show  that  FIM  I  and  FIM  0 
are  clearly  nonexponential  in  time,  while  FIM  X  is  close  to  an 
exponential.  A  stretched  exponential,  eq  19,  fits  FIM  I  better 
than  a  power  law,  eq  20,  at  pH  5.5  and  6.6.  Previously,  FIM  I 
was  fit  to  a  power  law  at  pH  7.^'  Because  the  deconvolution  of 
the  infrared  spectra  into  Ao,  A|,  and  A,  is  somewhat  ambiguous, 
the  details  of  the  time  dependence  are  difficult  to  pin  down.  To 
further  evaluate  the  data,  we  fit  the  temperature  dependence  of 
the  relaxation  functions  in  two  different  ways,  with  an  Arrhenius 
relation  and  with  eq  22.  The  results,  together  with  the  values  from 
ref  21,  are  given  in  Table  IV. 

Even  though  the  data  in  Table  IV  are  preliminary,  they  provide 
insight  into  the  relaxation  processes  in  Mb.  We  first  note  that 
FIM  X  differs  considerably  from  FIM  0  and  FIM  I.  Its 
preexponential  factor  as  obtained  from  an  Arrhenius  relation  is 


Figure  15  Relaxation  functions  4>(i.7T,  eq  18,  for  FIM  I,  FIM  X,  and 
FIM  0  FI.M  1  was  obtained  by  using  the  decrease  in  width  of  A).  FIM 
X  is  characterized  by  the  shift  of  the  peak  wavenumber  of  A«.  FIM  0 
is  characterized  by  the  change  in  the  area  of  Aj.  Sample:  swMbCO  in 
75%  glycerol/water  at  pH  6.6. 


TABLE  IV:  Parameters  Characterizing  the  Relaxation  Processes 
FIM  0,  FIM  I.  and  FIM  X* 


parameter 

pH 

FIM  0 

fast  slow 

FIM  1 

FI.M  X 

Time  Dependence 

d 

5.5 

6 

b 

0.3 

=  1 

6.6 

=  1 

0  2 

0.3 

=  1 

=  7 

=  1 

0.8 

c 

Temperature  Dependence-Arrhenius,  k 

=  A  exp(-£/R73 

log  (^/s  ') 

5.5 

b 

b 

20 

12 

66 

27 

70 

20 

13 

=7 

36 

£,  kj/mol 

5.5 

b 

b 

7S 

55 

6.6 

no 

280 

80 

55 

«7 

125 

Temperature 

Dependence,  k 

=  rlf  exp( 

Eq  22 

log  (.4b/s-') 

5.5 

b 

6 

8 

5 

6.6 

13 

34 

7 

5 

=7 

13 

21 

17 

£j,  kJ/mol 

5.5 

b 

b 

7 

6 

6.6 

9 

15 

7 

7 

=7 

10 

13 

9 

“The  time  dependence  of  the  relaxation  function  ♦(t.T)  described  by 
a  stretched  exponential,  eq  19,  with  parameter  d.  Temperature  de¬ 
pendence  given  in  terms  of  an  Arrhenius  (/4,£)  and  eq  22  (/4|,£|). 
FIM  0  has  a  fast  and  slow  component  for  pH  6.6  and  =»7.’'  FIM  X  is 
not  observed  for  pH  =  7.  Errors  are  ±25%.  ‘Parameters  not  well 
determined.  'FIM  1  for  pH  =»  7  is  better  described  by  a  power  law,  eq 
20  with  n  “  0.26. 

close  to  the  normally  expected  value  of  10'^  s"',  and  the  activation 
energy  could  place  it  into  tier  2.  Moreover,  it  is  close  to  exponential 
in  time.  We  therefore  believe  that  it  describes  a  less  complex 
phenomenon  than  either  FIM  0  or  FIM  1.  In  earlier  experiments^' 
at  a  higher  pH  we  did  not  observe  FIM  X;  it  may  consequently 
be  strongly  influenced  by  pH.  The  curves  in  Figure  13  hint  that 
more  detailed  studies  of  FIM  X  and  of  the  state  with  wavenumber 
1967.0  cm"'  are  feasible;  Near  180  K.  Mb  remains  in  this  state 
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for  a  sufnciently  long  time  so  that  additional  measurements  can 
be  performed. 

FIM  0  and  FIM  1  are  both  nonexponential  in  time.  FIM  I 
can  be  well  described  by  a  stretched  exponential.  FIM  0,  on  the 
other  hand,  is  more  complicated.  The  data  in  Figures  14  and  IS 
and  at  higher  pH^'  suggest  a  biphasic  behavior  in  which  at  least 
one  or,  possibly,  both  of  the  components  display  a  nonexponential 
time  dependence.  We  cannot  yet  decide  if  the  relaxation  is  due 
to  one  rather  complex  process  or  if  two  different  phenomena  are 
involved.  Moreo\er.  interference  from  FIM  X  is  possible.  To 
obtain  some  insight,  we  fit  <l>o  with  two  components  where  one 
is  a  stretched  exponential  while  the  other  is  close  to  exponential. 
These  components  are  labeled  “slow"  and  “fast"  in  Table  IV.  The 
temperature  dependencies  of  FIM  0  and  FIM  1  are  not  adequately 
described  by  an  Arrhenius  relation,  because  the  preexponentials 
take  on  values  of  10^  s"'  and  higher.  Similar  values  of  the 
preexponentials  arise  if  the  a  relaxation  in  glasses  is  characterized 
by  an  Arrhenius  relation.  In  glasses,  where  the  a  relaxation  can 
be  observed  over  a  wide  range  in  T,  it  is  known  that  an  Arrhenius 
plot  does  not  yield  a  straight  line  and  that,  for  instance,  eq  21  or 
22  must  be  used.  We  invoke  the  analogy  to  glasses  and  use  eq 
22  to  parametrize  the  data  in  Table  IV.  The  values  of  the 
preexponential  and  energy  for  FIM  1  appear  to  depend  strongly 
on  pH.  We  assume,  again  in  analogy  to  glasses,  that  both  FIM 
0  and  FIM  1  describe  large-scale  motions  in  the  protein.  A 
complete  characterization  of  relaxation  in  MbCO  is  dimcult  owing 
to  the  fact  that  several  relaxation  processes  occur  in  a  narrow- 
temperature  interval  around  180  K.  In  addition,  the  biphasic 
nature  of  FIM  0  is  not  understood.  Therefore,  all  relaxation 
processes  in  MbCO  will  have  to  be  studied  in  more  detail  before 
they  can  be  fully  interpreted. 

8.  Pressure  and  Protein  Reactions 

We  discuss  here  one  simple  example,  the  binding  of  CO  and 
Oj  to  myoglobins  at  low  temperatures,  but  believe  that  the  general 
concepts  will  be  needed  to  explain  the  high-temperature  data  and 
are  also  applicable  to  other  reactions  in  complex  systems. 

8.1 .  Pressure  Effect  on  Ligand  Binding — A  Puzzle'i  The 
measurements  of  the  effects  of  pressure  on  the  association  rates 
of  CO  and  Oj  to  Mb  and  hemoglobin  (Hb)  by  HasinofP*  gave 
a  puzzling  result.  The  activation  volumes  [eq  4]  for  the  binding 
of  O2  were  positive  (+8  cmVmol  for  Mb,  +5  cmVmoI  for  Hb) 
while  those  for  CO  were  negative  (-9  cm’/mol  for  Mb,  -3  and 
-21  cmVmol  for  the  “fast"  and  “slow"  reaction  in  Hb).  The 
negative  values  observed  for  CO  conform  to  expectation:  The 
binding  of  CO  involves  formation  of  a  covalent  bond  between  CO 
and  the  heme  iron,  and  such  bond  formation  is  typically  char¬ 
acterized  by  an  activation  volume  of  about  -10  cmVmol.”  In 
contrast,  the  positive  activation  volume  for  Oj  is  surprising.  The 
puzzle  is  magnified  on  closer  inspection.  The  binding  of  CO  and 
Oj  is  controlled  at  the  heme  iron,  and  this  binding  step  is  similar 
for  both  ligands.*®  Moreover,  binding  of  CO  to  heme  has  an 
activation  volume  of  +2  cm’/mol.*'  These  values  demonstrate 
that  the  binding  process  is  more  complex  than  a  simple  one-step 
reaction  and  that  the  protein  must  play  a  role  in  the  control  of 
the  binding  reaction.  We  will  show  here  that  the  effect  of  pressure 
on  the  simple  association  reaction  is  sufficiently  complex  that 
experiments  performed  only  at  physiological  temperatures  and 
using  only  observation  in  the  Soret  are  incapable  of  elucidating 
the  mechanisms  of  such  biomolecular  reactions. 

8.2.  Low-Temperature  Binding  of  CO  and  O2  to  Mb.  In  order 
to  discuss  the  effect  of  pressure  on  the  binding  of  ligands  to  heme 
proteins,  the  general  features  of  the  binding  process  must  be 
described  first.  While  pressure  data  over  the  entire  temperature 
range  from  4  to  300  K  exist,*^  we  restrict  the  discussion  here  to 
the  low-temperature  regime  below  about  160  K.  This  restriction 

(58)  Hasinofr,  B.  B.  Biochemistry  1974,  13,  3111. 

(59)  Van  Eldick.  R.;  Asano,  T.:  U  Noble.  W  J  Chem  Rev  19B9, 89.  549. 

(60)  Frauenfclder.  H.;  Wolynes,  P  Science  1985.  2!9,  337. 

(61 )  Caldin,  E  F.;  Haainoff,  B.  B.  J.  Chem.  Soc.,  Faraday  Trans  I  1975. 
71.  515 

(62)  Alberding.  N.  A.  Ph.D  Thesis.  University  of  Illinois  at  Urbana- 
Champaign,  I97g. 
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Figure  16  Rebinding  of  CO  (a)  and  O;  (b)  to  sperm  whale  myoglobin 
after  flash  photolysis.  AA(t).  the  absorption  change  at  time  1  after  the 
laser  flash,  is  proportional  to  the  fraction  of  Mb  molecules  that  have  not 
rebound  a  ligand  The  data  have  been  scaled  to  A.4(0)  for  swMbCO  in 
the  F  state  Solvent:  75%  glycerol/water,  pH  7,0  The  pressure  states 
F,  FP.  PF,  and  PFR  are  defined  in  Figure  8 

permits  us  to  concenuan?  or.  iiie  central  aspects  of  the  binding 
reaction 

The  binding  of  CO  and  Oj  to  heme  proteins  at  temperatures 
below  about  160  K  has  been  studied  in  detail  by  using  flash 
photolysis. The  sample,  for  instance  MbCO,  is  brought  to 
the  desired  temperature.  The  bond  between  the  CO  and  the  heme 
iron  is  broken  by  a  laser  flash.  The  CO  moves  into  the  heme 
pocket  and  rebinds  from  there.  We  denote  the  bound  slate  as  A 
and  the  state  with  the  ligand  in  the  heme  pocket  by  B  (Figure 
lb).  Rebinding  B  -•  A  is  monitored  by  observing  the  corre¬ 
sponding  spectral  changes,  either  in  the  Soret'*  or  the  infrared." 
We  denote  with  N(t)  the  survival  probability,  the  fraction  of  Mb 
molecules  that  has  not  rebound  a  ligand  at  the  time  t  after 
photodissociation.  N(t)  is  nonexponential  in  time  and  can  usually 
be  approximated  by  a  power  law,  eq  20.'*  Examples  of  such 
nonexponential  rebinding  are  given  in  Figures  16  and  17. 

The  nonexponential  rebinding  is  explained  by  postulating  that 
each  protein  molecule  is  in  a  different  conformational  substate 
and  that  the  different  substates  possess  different  activation  en¬ 
thalpies  for  the  binding  transition  B  -*  A.'*  To  describe  the 
observed  N(i),  we  write 

N(t)  =  f  dHg„  g(//gj  expI-A-fWgJf]  (23) 

Here  g{Hgt,)  dZ/gA  is  the  probability  of  finding  a  protein  molecule 
with  an  enthalpy  barrier  between  Wra  and  +  dWgA.  and 
*(WgA)  is  the  rate  coefficient 

lr(//gA)  =  zlgA  espl-Hg/^/ RT]  (24) 

The  preexponenlial,  /4gA.  is  usually  taken  as  /4gA  =  A{T/To), 
where  the  reference  temperature  To  is  100  K  and  A  is  temperature 
independent. 

8.3.  Pressure  Effects  Observed  in  the  Soret  Band.*^  We  expect 
pressure  to  have  two  dominant  effects  on  reaction  rates,  similar 

(63)  Ansari,  A.,  Di  lorio,  E  E.:  DIoll.  D  D.;  Frauenfeldcr,  H.;  Iben.  1 
E  T..  Langcr.  P  ;  Roder,  H  :  Saukc,  T  B  ;  Shyamsunder.  E  Biochemistry 
1986,  25.  3139 
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Figure  17.  Rebinding  of  CO  (a)  and  Oj  (b)  to  horie  mvogl  'bin  afier 
pholodissocialion.  The  data  have  been  scaled  to  A/1(0)  for  swMbCO  in 
the  F  state.  Solvent  TST  glycerol/water,  pH  7  0 

to  the  two  effects  on  spectral  line  shifts.  In  each  given  confor¬ 
mational  substate,  the  reaction  rate  will  be  changed  through  the 
activation  volume  according  to  eq  3  and  4.  Pressure  can.  however, 
also  shift  the  populations  among  substates  and  thus  affect  the 
reaction  rate.  To  separate  the  two  phenomena,  we  use  the 
technique  shown  in  Figure  8  and  described  in  section  6  for  spectral 
shifts. 

The  pressure  experiment  is,  in  principle,  straightforward.  W'e 
observe  rebinding  after  flash  photolysis  for  the  four  pres.sure  states 
F.  FP,  PF,  and  PFR  of  Figure  8  for  sperm  whale  (sw  Mb)  and 
horse  (hMb)  myoglobin.  In  practice,  low-temperature  high- 
pressure  flash  photolysis  experiments  lead  to  many  experimental 
frustrations,  and  the  help  and  guidance  of  Harry  Drickamer  has 
been  essential  for  this  work.  A  selection  of  experimental  results 
are  shown  in  Figure  16  and  1 7,  where  log  JsA(t)  is  plotted  versus 
log  i.  The  absorbance  change  SA(t)  at  time  t  is  proportional  to 
the  survival  probability  N{t).  eq  23.  The  binding  of  CO  was 
monitored  at  436  nm  and  that  of  O2  at  440  nm. 

Even  without  detailed  evaluation,  the  data  in  Figures  16  and 
1 7  show  that  both  activation  and  conformation  effects  are  pre,sent, 
that  sperm  whale  and  horse  Mb  behave  differently,  and  that  CO 
and  Oj  also  display  clear  differences. 

8.4.  Evaluation  of  the  Pressure  Data.  Evaluation  entails  three 
aspects;  the  nonexponential  time  dependence,  extraction  of  ac¬ 
tivation  volumes,  and  interpretation. 

The  nonexponential  time  dependence  is  described  in  terms  of 
the  probability  density  giHef)  in  eq  23.  In  principle,  we  can  invert 
eq  23  for  each  pressure  state  and  find  the  distributions  for  the 
four  pressure  states.^  However,  this  procedure  does  not  yield 
activation  volumes  directly  and  is  technically  often  difficult  be¬ 
cause  IV(/)  cannot  easily  be  measured  over  a  sufficiently  extended 
time  range.  We  introduce  here  a  simple  approach  that  gives  some 
of  the  essential  information  directly.  Consider  the  nonexponential 
rebinding  curve  and  the  corresponding  probability  density  dis¬ 
tribution  shown  in  Figure  18.  Equation  23  shows  that  at  a  given 
time  r  only  a  small  band  of  activation  enthalpies  contributes  to 
the  rebinding.  Proteins  with  rates  satisfying  k(Hg^)T  »  I  have 
already  rebound  their  ligand  while  those  with  k(Hg^)T  «  I  rebind 
much  later.  Only  proteins  with  barriers  centered  at 

- /?rin  M,*rl  (25) 

contribute  to  binding  at  time  r.  We  use  a  temper.-iture-independent 
preexponential  e<l  25  for  simplicity.  Equation  25  assigns 
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Figure  18  S'oncxponential  rebinding.  N(l).  and  ihe  corresponding 
probdbiliiy  density  distribution.  gtHg*).  At  the  time  r  only  proteins  with 
enthalpies  near  rebind. 

an  average  activation  enthalpy  Wba(’’)  I'd  «dch  time  r. 

This  approach  permits  us  to  find  an  approximate  value  for  the 
activation  volume:  Pressure,  as  indicated  in  Figure  18,  can  change 
the  rebinding  If  no  conformational  changes  occur,  the  shift  from 
r  to  r'  at  constant  ,V(  r)  is  caused  by  an  activation  volume.  The 
pressure  dependence  of  the  rate  coefficient  is  given  by  eq  3  and 
4  as  k{P)  =  const  X  cxp[-FF*8^//f71.  If  k{P)  is  measured  at 
two  pressures.  P  and  P'.  the  activation  volume  is  obtained  from 
the  ratio  k(P)/k(P^  as 

D  J" 

f"BA  =  (26) 

With  k{P)  =  l/r(P),cq  26  permits  the  immediate  evaluation  of 
the  activation  volume  from  the  data  in  Figures  16  and  17. 
Moreover,  by  obtaining  F^ba  ^t  various  values  of  N(t),  we  find 
the  dependence  of  the  activation  volume  on  the  barrier  enthalpy 
at  the  initial  pressure.  Since  the  barrier  enthalpy  characterizes 
the  substates,  we  are  able  to  investigate  the  substate  dependence 
of  k'*BA 

If  pressure  causes  a  conformational  change,  the  situation  is 
different.  The  two  curves  labeled  P  and  P'  in  Figure  18  then  can 
correspond  to  the  same  external  pressure,  but  to  different  con¬ 
formations  The  curves  labeled  F  and  PFR  in  Figures  16  and  17 
represent  such  a  pair.  We  can,  nevertheless,  use  eq  26  to  define 
a  conformation  volume  to  describe  the  rate  change  in  going 
at  a  fixed  value  of  N(t)  from  r  to  P.  The  pressures  Pand  F'then 
do  not  correspond  to  the  external  pressure,  but  to  the  pressures 
at  which  the  sample  was  frozen.  Rebinding  data  taken  at  the  four 
points  F,  FP,  PF,  and  PFR  in  the  PT  plane  (Figure  8)  thus  permit 
the  definition  of  three  activation  volumes 

P^l  =  ^FP-F-  *''*2  =  f^PF-PFR-  =  ^PFR-F  (27) 

The  subscripts  indicate  between  which  pairs  of  curves  the 
difference  tin  k(P)  -  In  k(P')\  has  to  be  taken.  The  pressure 
difference  for  all  volumes  is  the  pressure  at  the  points  0  and  P 
in  Figure  8.  F^i  and  F*2  arc  activation  volumes.  Vq,  the  con¬ 
formation  volume,  is  a  parameter  with  units  of  volume,  but  with 
no  direct  interpretation  in  terms  of  a  real  volume.  Vq  expresses 
the  direction  and  effect  of  conformational  changes. 

Values  of  the  activation  and  conformation  volumes,  obtained 
with  eq  26  from  the  data  in  Figures  16  and  17,  are  given  in  Table 
V.  N(t)  is  defined  in  Figure  18.  is  obtained  from  eq 

25.  with  temperature-independent  preexponentials  obtained  by 
using  a  gamma  distribution  for 

The  values  in  Table  V  show  that  horse  and  sperm  whale 
myoglobin  behave  very  differently  under  pressure  even  though 
the  rebinding  curves  at  F  =  0  look  very  similar.  In  both  proteins, 
the  activation  volumes  for  MbCO  arc  negative  and  of  the  expected 

(64)  Young.  R.  D .  Bowne.  S.  F  J  Chem.  Phyi.  I9S4,  81.  3730. 
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TABLE  V:  Actiittion  Volumes  V*,  and  V*]  and  Conformation 
Volume  Describin):  the  Effect  of  Pressure  on  the  Four  States  of 
Fifture  8* 


sample 

T.  K 

A.4(t) 

(OD) 

kJ/mol 
(state  F) 

1^1. 

cm’/mol 

cm'/mol 

V’c. 

cm’/ mol 

swMbCO 

80 

03 

12 

-8 

-7 

-8 

0.1 

16 

-7 

-7 

-9 

120 

0  1 

16 

-9 

-7 

-9 

0.01 

21 

-6 

-7 

-9 

swMbOj 

80 

0  3 

12 

-2 

0 

+7 

120 

0.3 

12 

-2 

-7 

+7 

0.1 

16 

0 

-3 

+  13 

hMbCO 

80 

0.3 

13 

-5 

-10 

+4 

120 

03 

12 

-6 

-10 

+5 

0.1 

16 

-5 

-8 

+3 

0.01 

21 

-3 

-6 

+2 

hMbOj 

80 

0  3 

13 

0 

-4 

-8 

120 

0  3 

13 

+  2 

-1 

-7 

0  1 

18 

-1-4 

+4 

-7 

•  Temperature- 

independent  preexponentials. 

Aba:  1-2 

X  10’  s-' 

(swMbCO),  1.8  X  10’ s-'  (swMbO;),  2  0  x  10’  s  '  (hMbCO),  1.8  x 
10’  s''  (hMbOj).  Rebinding  monitored  in  Sorel  band;  MbCO  at  436 
nm,  MbOj  at  440  nm  Samples  sperm  whale  (swMb)  and  horse 
(hMb)  myoglobin  in  75*^  glycerol/water  (v/v)  with  phosphate  buffer 
pH  7.0.  Error  in  F  is  ±1  cm’/mol. 

Table  VL  Activation  and  Conformation  Volumes  of  the  Substates 
Kf  and  Aj  for  Binding  of  CO  to  Sperm  Whale  Myoglobin'' 


substale 

pH 

r.  K 

A(r) 

kJ/mol 
(state  F) 

em'/mol 

cm’/mol 

cm’/mol 

sum 

5.5 

25 

0.5 

7 

-4 

-16 

Aq 

0  6 

7 

0 

-12 

sum 

66 

40 

0.5 

11 

-10 

-2 

-19 

Aq 

0  3 

1 1 

-5 

-2 

-11 

A. 

06 

10 

-8 

-2 

-5 

sum 

66 

50 

0.3 

13 

-8 

-3 

-15 

Aq 

0  1 

15 

-4 

-3 

-11 

A, 

0  3 

13 

-9 

-s 

-s 

*  Preenponentials  for  Aoand  A,  from  ref  II  The  50  K  data  arc  not 
shown  in  the  figures  Error  in  Pis  ±1  cm’/mol  Solvent:  7571  gly¬ 
cerol/water  (v/v)  and  potassium  phosphate  buffer. 

magnitude  for  bond  formation;  pressure  speeds  up  CO  rebinding 
The  magnitude  of  the  activation  volume  for  CO  decreases  with 
increasing  barrier  height.  The  activation  volume  for  rebinding 
of  Oj  is  small  for  both  swMb  and  hMb,  but  clearly  negative  for 
swMb  and  mostly  positive  for  hMb  These  values  for  dioxygen 
solve  part  of  the  puzzle  created  by  the  room  temperature  mea¬ 
surements:  The  activation  volume  for  the  rebinding  step  at  the 
heme  is  not  large  and  positive,  but  small  and  mainly  negative.  The 
conformational  effects,  which  are  characterized  by  Pc-  ^me  large 
and  negative  for  swMbCO  and  hMbOj,  large  and  positive  for 
swMbOj,  and  small  and  positive  for  hMbCO.  The  large  difference 
between  the  conformation  volumes  for  swMbCO  (-9  cm’/mol) 
and  swMbOj  (+7  to  +13  cm’/mol)  obtained  from  the  Sorel  band 
suggests  a  strong  interplay  between  ligand  and  protein. 

S.5.  CO  Rebinding  Observed  in  the  Infrared  Spectrum.  The 
Soret  data  show  that  protein  structure  and  ligand  properties  affect 
the  pressure  dependence  of  ligand  binding,  but  they  cannot  tell 
us  which  tier  of  substates  is  responsible.  This  information  is 
obtained  by  monitoring  the  ligand  rebinding  using  the  infrared 
CO  stretch  bands.  Rebinding  data  after  flash  photolysis  for  the 
bands  Aq  and  A)  are  given  in  Figures  19  and  20.  At  pH  6.6, 
both  bands  can  be  measured:  at  pH  5.5,  the  A,  band  in  the  states 
PF  and  PFR  is  loo  small  to  be  measured.  Figures  19  and  20  show 
that  pressure  affects  the  rebinding  in  both  substates  but  that  the 
changes  differ  considerably  for  the  two  bands  From  the  measured 
curves,  the  activation  volumes  P*,  and  and  the  conformation 
volume  Pc  are  determined  by  the  method  described  in  subsection 
8.4;  the  results  are  summarized  in  Table  VI.  The  data  in  Table 
VI  and  the  curves  in  Figure  19  show  that  the  speed-up  of  the  CO 
rebinding  has  two  causes:  Pressure  shifts  the  population  from  the 


log  ( l/i) 

Figure  19  Rebinding  of  CO  to  swMb  after  photodissociation  at  40  K 
as  a  function  of  time  for  different  pathways  (see  Figure  8).  (a)  Re¬ 
binding  of  the  total  A  substates;  (b)  rebinding  to  A®:  (c)  rebinding  to  A,. 
Sample.  swMbCO  in  glycerol/water  (v/v)  at  pH  6.6. 


Figure  20  Rebinding  of  CO  to  swMb  after  photodissociation  at  25  K 
as  a  function  of  lime  for  different  pathways  (see  Figure  8).  (a)  Re¬ 
binding  of  the  total  A  subslales;  (b)  rebinding  to  A^  Sample  swMbCO 
in  75^  glycerol/water  (v/v)  at  pH  5.5. 

slower  A,  substate  to  the  faster  Ag  subslate,  and  the  rebinding 
of  both  subslates  speeds  up.  The  rebinding  of  each  substate  is 
nonexponential  in  time.  These  results  require  a  generalization 
of  eq  23  to  include  multiple  A  substates: 

m  =  LH'ifP)  J dHsa  «,(«ba)  cxp[-k,{H^^)t]  (28) 

The  sum  extends  over  the  substates  of  tier  0,  w.fP)  is  the  weight 
of  substate  A,  at  pressure  P,  and  is  the  corresponding 

probability  distribution  of  the  barrier  height.  The  rale  coefficients 
/i,(Hb^)  are  given  by  eq  24,  but  the  preexponential  Ag*  may  be 
different  for  each  substate."  The  speed  up  for  CO  rebinding  to 
swMB  is  described  by  the  changes  in  the  weights  k’,(P)  and  in 
the  distributions 

The  two  substates  Ag  and  A,  behave  differently  under  pressure 
as  indicated  by  the  parameters  in  Table  VI.  The  increased  re- 
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Figure  21  Rebinding  of  CO  to  swMb  after  phoiodissociaiion  at  40  K 
as  a  function  of  time  at  different  wavenumbers.  (a|  State  F.  (b)  state 
FP  Sample:  swMbCO  in  75'’c  glycerol/water  (v/v)  at  pH  6  6.  (O) 
1941  cm  ',  (O)  area  of  A,,  (□)  1951  cm"' 


and  B|  -•  A,,  can  occur  **  At  high  temperatures,  the  two  types 
of  transitions  may  have  about  equal  rates.  Here  we  discuss  low 
temperatures,  T  <  T",,,  where  conformational  transitions  are 
absent,  as  was  discussed  in  sections  5  and  7.  The  survival  prob¬ 
ability  iV(r)  of  the  photodissociated  species  can  then  be  written 
as 

,V(f)  =  wo(/>)  expl-M/’)'l  +  *^i(n  enpl-^iWtl  (29) 

Here  wq{P)  and  wilP)  are  the  substate  populations,  and  kQ(P) 
and  k,{P)  are  the  rate  coefficients  for  CO  binding  from  subslates 
Bo  and  B|.  The  survival  probabilities  .V(r)  at  time  r  (Figure  18) 
for  the  four  pressure  states  of  Figure  8  are 

Nfir)  =  H'o(O)  exp|-<:o(0)T|  +  w-iiO)  expl-<i|(0)T| 

Apytr)  =  wo(0)  exp|-l:o(/’)T|  +  >V|(0)  exp|-/C|(/’)r| 

.Vpy(r)  =  wo(/’)  exp|-lro(^)’'l  +  *'i(^)  expl-<:|(P)r| 

■''pFR(r)  =  Wolf)  >:>‘Pl'^o(())’’l  +  "’^P)  expl-<:|(0)r|  (30) 

Following  the  arguments  in  section  6.1,  we  introduce  three  in¬ 
dependent  differences  With  the  normaliz.ition  y^olP)  +  M;(P) 
=  I,  and  the  abbreviations 

6,(r)  =  exp[-<:  (Pfr]  -  exp[-A:,(0)rl,  /  =  0,  I  (31) 


binding  rate  of  .Aq  is  caused  by  a  large  conformation  volume,  while 
the  increase  in  A,  is  the  result  of  a  large  activation  volume 
Pressure  thus  mainly  affects  go(^BA)  4nd  The  difference 

between  .Aq  and  A,  is  also  evident  when  the  rebinding  is  monitored 
at  different  wavenumbers  within  the  same  A  band.  Each  substatc 
of  tier  0  is  comprised  of  a  large  number  of  substates  of  tier  1  and 
is  inhomogeneously  broadened.  If  these  CS'  have  different  ac¬ 
tivation  energies,  different  wavenumbers  within  the  same  band 
should  rebind  differently  at  P  =  0  (state  F).  This  effect  is  called 
kinetic  hole  burning.'*’  ” ”  Figure  21a  shows  that  this  prediction 
is  indeed  true  for  A,.  If  different  CS'  also  have  different  activation 
volumes,  rebinding  in  the  state  FP  should  differ  from  that  in  F. 
Figure  21b  shows  that  the  rate  increase  from  1941  to  1951  cm  ' 
is  indeed  larger  in  FP  than  in  F.  In  substate  A,,  both  activation 
and  conformation  effects  are  important.  In  contrast,  kinetic  hole 
burning  is  very  small  in  Aq. 

In  Table  VI,  the  activation  and  conformation  volumes  are  listed 
for  the  sum  of  the  A  substates.  We  v.ould  expect  the  sum  of  the 
infrared  bands  to  give  the  same  parameters  as  the  Soret  band 
While  the  two  sets  show  the  same  trend,  there  are  quantitative 
differences.  V*2-  in  particular,  is  much  smaller  in  the  infrared 
than  in  the  Soret  band.  The  different  measurement  techniques 
may  be  responsible  for  the  discrepancy;  In  the  IR  regions  the 
band  areas  are  monitored;  in  the  Soret  region,  only  a  small 
wavelength  band  is  observed.  The  discrepancy  thus  leads  to  a 
prediction,  not  yet  tested:  The  pressure  effect  observed  should 
depend  on  wavelength  in  the  Soret  region,  just  as  in  the  IR  region. 

The  MbCO  results  obtained  in  the  infrared  region  permit  some 
indirect  conclusions  about  the  binding  of  dioxygen  The  stretch 
bands  of  Oj  are  difficult  to  observe  and  make  the  IR  monitoring 
of  Oj  binding  nearly  impossible.  However,  the  data  in  Table  V 
show  that  Vq  fot  Oj  binding  to  swMb  is  large  and  positive; 
Pressure  slows  down  Oj  binding  and  the  main  cause  for  the 
slowdown  is  a  conformational  shift  ar.d  not  an  activation  volume. 
In  contrast,  the  binding  of  Oj  to  horse  myoglobin  is  more  weakly 
affected  by  pressure,  and  the  conformation  volume  has  the  opposite 
sign  from  swMb,  It  is  amusing  to  speculate  whether  the  difference 
is  related  to  the  fact  that  whales  dive  and  horses  do  not. 

S  6.  A  Model  of  Acliuation  and  Conformation  Effects.  The 
evaluation  introduced  in  subsection  8.4  is  only  semiquantitative. 
A  full  treatment  of  the  data  is  very  complicated,  and  we  will  not 
perform  it  here.  We  will,  however,  sketch  the  approach  for  the 
simple  situation  shown  in  Figure  1.  We  assume  two  conforma¬ 
tional  substates  of  tier  0,  labeled  0  and  I .  The  protein  then  has 
four  states,  the  bound  substate  Aq  and  A|  and  the  two  photo- 
dissociated  pocket  states  Bo  and  B,.  The  notation  for  the  Bq 
substate  differs  from  previous  usage."  Two  conformational 
transitions,  Ao  A,  and  Bo  B,.  and  two  reactions.  Bo  Ao 


we  obtain 

A,'(r)  =  .V,p(T)  -  .\p(r)  =  K,,(0)  3o(r)  -F  k,(0)  (5,(r) 

Aj'tr)  =  .Ap|.(r)  -  .Vp|.p(r)  =  vvqIP)  '^oI  r)  +  ^\\P)  5|(r) 
A^(r)  =  .Vpiplr)  -  ,\f(r)  = 

(wiiP)  K|(nil[exp|-k|(0)r|  ■  exp|-k(,(0)Tl]  (32) 
As  in  cq  14.  we  also  form  the  difference 

A;'(r)  -  A,'(r)  =  |>V|(P)  -  H|(0)1|6;I  r)  -  5otr)l  (33) 

Similar  arguments  to  those  in  section  6,1  show  the  following 
nonvanishing  A,'  and  A.'  prove  that  the  activation  volumes  for 
the  substates  are  nonzero,  nonvanishing  A/  shows  that  pressure 
changes  the  population  of  the  substates  and  that  the  rates  in  the 
two  subslates  are  dil'fcrcni;  and  a  nonvanishing  difference  Aj'  - 
A|'  proves  that  the  activation  volumes  affect  6,(r)  for  the  two 
substates  differently. 

While  the  approach  given  so  far  describes  some  of  the  main 
features  of  the  observed  pressure  effects,  it  must  be  generalized 
in  at  least  two  ways:  cq  29  assumes  that  rebinding  is  exponential 
in  time  for  each  substatc  To  'ake  the  nonexponential  rcbtnding 
into  account,  cq  28  must  be  used.  The  second  generalization  is 
demanded  by  the  hierarchy  of  substates.  In  eq  30,  the  same  rate 
coefficients  k^iP)  and  k^(P)  are  used  for  the  two  pressure  states, 
FP  and  PF.  The  transitions  Bp  -*  Ap  and  B,  -•  A,  involve, 
however,  a  large  number  of  transitions  between  substates  of  tier 
I.  Pressure  changes  the  distribution  of  these  CS'.  and  the  rates 
for  each  pressure  state  (FP.  PF,  F.  and  PFR)  can  be  different. 
This  difference  is  evident  in  Table  VI  where  the  parameters  F"*, 
and  F’*2  arc  not  identical.  The.se  generalizations  exceed  the  scope 
of  the  present  paper,  and  we  consequently  do  not  complete  the 
evaluation  here. 

9.  Summary  and  Outlook 

As  a  thermodynamic  variable,  pressure  is  as  important  as 
temperature  in  the  study  of  proteins.  In  the  present  paper,  we 
have  treated  five  areas  where  pressure  is  an  essential  tool  (section 
1 ).  In  each  of  the  five  areas,  the  pressure  experiments  have  yielded 
new  insights  and  opened  new  avenues.  Pressure  affects  proteins 
differently  than  it  does  crystals  because  proteins  are  complex 
systems  and  can  assume  a  large  number  of  different  conforma¬ 
tional  substates,  arranged  in  a  hierarchy  of  tiers  (section  2) 
Pressure  consequently  affects  proteins  in  two  ways;  it  can  act  on 
a  protein  in  a  given  substate  (elastic  effect),  and  it  can  shift  the 
distribution  of  substates,  thereby  influencing  protein  structure, 
dynamics,  and  function  (conformational  effects). 
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Myoglobin  shows  a  glass  transition  near  180  K  if  embedded 
in  a  glycerol/water  solvent.  Well  above  the  transition  temperature, 
protein  substates  are  in  equilibrium  and  the  thermodynamic  pa¬ 
rameters  describing  substates  of  tier  0,  in  particular  their  volume 
differences,  can  be  determined  (section  4).  Our  results  imply  that 
the  internal  energy,  the  entropy,  and  the  volume  of  myoglobin 
are  functions  of  temperature.  The  experiments  also  show  that 
infrared  measurements  over  a  wide  range  in  temperature  are 
essential  for  protein  studies;  if  data  are  taken  only  in  the  visible 
region  or  over  a  small  temperature  range,  an  unambiguous  in¬ 
terpretation  is  difficult  or  impossible. 

Below  about  180  K,  myoglobin  becomes  glasslike,  the  me¬ 
tastability  is  proven  by  combining  pressure  and  temperature 
variations  and  demonstrating  that  the  infrared  spectrum  of 
myoglobin  depends  on  the  path  in  the  TP  plane  by  which  a  final 
temperature-pressure  point  is  reached  (section  5) 

The  technique  of  taking  different  pathways  to  a  given  TP  point 
below  the  glass  temperature  also  permits  the  separation  of  elastic 
and  conformational  pressure  effects  on  spectral  bands  (section 
6)  These  experiments  show  that  different  substates  of  tier  0 
possess  different  elastic  and  conformational  properties 

One  of  the  goals  of  protein  studies  is  the  measurement  and 
description  of  protein  motions  The  pressure  jump  technique 
provides  a  generally  applicable  approach  to  explore  protein  motions 
(section  7);  it  shows  that  relaxation  phenomena  are  complicated 
but  can  be  unravelled 

The  most  impi'rtant  use  of  pressure  may  well  be  for  the  study 
of  protein  reactions  Here  again  the  different  pathways  in  the 
PT  plane  permit  the  separation  of  the  effects  of  pressure  into  an 
activation  and  a  conformation  component,  characterired  b\  ac¬ 
tivation  and  conformation  volumes  (section  8)  A  surprising  result 
of  these  studies  is  the  large  difference  m  behavior  between  closely 
related  proteins  such  as  sperm  whale  and  horse  myoglobin  \\  hile 
It  IS  dangerous  to  extrapolate  from  a  small  sample,  it  appears  that 
the  pressure  parameters  give  as  much  or  more  information  than 
activation  enthalpies 

The  experiments  and  models  treated  in  the  present  paper  are 
only  a  beginning  and  call  for  a  large  number  of  additional  ex¬ 
periments  and  for  better  theoretical  treatments  Foremost  is  the 
extension  of  the  experiments  to  other  proteins,  both  wild-type  and 
sy  nthetic  mutants  VN'ithout  a  larger  data  base  the  results  cannot 
be  fully  interpreted  It  is  not  clear  which  phenomena  arc  general 
and  which  ones  are  specific  (or  accidental)  in  a  particular  protein 
F.qually  important  is  the  connection  of  the  pressure  data  to 
structural  and  calorimetric  measurements  The  experiments 
described  here  show  that  the  dynamic  properties  of  myoglobin 
depend  strongly  on  pH  and  suggest  that  even  the  overall  structure 
and  the  specific  heat  may  differ  in  the  different  substates  of  tier 
0  X-ray  and  specific  heat  data  as  a  function  of  pH  can  verify 
or  discredit  this  mode! 

An  ordinary  glass  usually  shows  one  glass  transition  (o  re¬ 
laxation)  with  non-.Arrhenius  temperature  dependence,  secondary 
glass  transitions  (d.y)  are  less  pronounced  Myoglobin  shows 
at  least  three  relaxation  phenomena  near  180  K  Two  have 
properties  characteristic  of  a  glass  transition,  and  the  other  may 
be  similar  to  the  8  relaxation  Such  a  behavior  is  not  entirely 
unexpected:  Different  domains'*  in  the  protein  could  well  have 
related  but  somewhat  different  relaxation  phenomena  The 
technique  discussed  in  section  7  will  permit  more  studies  of  these 
relaxation  processes  by  looking  not  only  at  the  CO  stretch  bands, 
but  also  at  other  IR  bands 

The  sensitivity  of  the  rebinding  to  swMb  and  hMb  to  the 
activation  and  conformation  volumes  suggests  another  class  of 
experiments;  Monitoring  rebinding  in  the  IR  and  the  use  of 
different  tsotopes  of  C  and  O  earlier  permitted  us  to  conclude  that 
CO  does  not  bind  like  a  point  particle,  but  must  undergo  some 
rotational  motion  The  combination  of  isotopic  replacements 
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with  rebinding  measurements  at  various  points  in  the  TP  plane 
should  provide  an  even  closer  look  at  the  binding  step  at  the  heme 
iron 

in  the  present  paper  we  discussed  only  the  binding  of  CO  and 
O2  below  180  K,  where  the  entire  binding  process  takes  place  inside 
the  heme  pocket  We  have  performed  measurements  of  the 
pressure  dependence  of  binding  up  to  about  300  K  These  data 
again  provide  additional  information  on  the  rather  complicated 
binding  process,  b’lt  discussion  of  the  results  is  beyond  the  scope 
of  this  paper. 

All  experiments  treated  in  the  present  paper  share  one  problem 
(and  opportunity)  not  yet  discussed  in  detail  The  data  obtained 
all  refer  not  to  the  properties  of  an  isolated  protein  but  to  the 
protein-hydration  shell-solvent  system  This  observation  is 
particularly  important  in  pressure  experiments,  because  pressure 
changes  the  viscosity,  the  relaxation,^  the  pH,  and  the  pA'  of  the 
solvent  The  results  of  the  present  work  show  that  the  properties 
of  Mb  also  depend  strongly  on  these  parameters  It  will  therefore 
be  important  ti)  study  the  crucial  features  of  the  energy  landscape, 
the  protein  dynamics,  and  the  protein  reactions  with  different 
solvents  in  order  to  understand  these  effects  As  one  example  we 
mention  that  we  have  evaluated  the  binding  data  for  Mb  with  the 
Arrhenius  relation,  eq  24.  with  the  pree’ponenlial  factor  .-Ig, 
assumed  to  be  pressure  inoependent  Kramers’  theory  implle.^. 
however,  that  the  prefaclor  depends  on  viscosity, “  and  viscosity 
IS  prcs^ure  dependent  Pressure  thus  affects  the  reaction  rate  not 
only  through  activation  and  conformation  volume  but  also  through 
changes  in  the  precxponential  factor,’'’  and  this  change  should 
also  be  studied  It  is  not  unfortunate  that  we  do  not  obtain 
information  about  the  isolated  protein  because  the  functioning 
protein  in  a  living  system  is  not  in  a  standard  aqueous  solution 
but  IS  alwavs  part  of  a  more  complex  system  Studies  with  dif¬ 
ferent  solvents  and  .solvent  conditions,  therefore,  can  tell  us  much 
about  possible  effects  of  the  environment  on  proteins  within  living 
systems 

ThcorcUcal  work  is  indicated  in  at  least  two  different  directions 
In  section  8,  we  introduced  a  highly  simplified  treatment  of  the 
effect  of  pressure  on  the  reaction  rate  This  treatment  must  be 
extended  to  take  the  nonexpsincntial  time  dependencies  and  the 
hierarchy  of  substates  into  account  A  more  difficult  theoretical 
problem  ts  the  connection  between  the  binding  process  and  mo¬ 
lecular  dyn.imics  com, putations  '  L  Itimately,  this  connection 
may  elucidate  the  way  in  which  protein  structure  and  dynamics 
control  protein  function 

Pressure  effects,  to  which  Harry  Drickamcr  has  devoted  es¬ 
sentially  his  entire  .scientific  work,  promise  to  be  as  important  and 
fruitful  in  the  study  of  proteins  as  they  have  been,  and  still  are. 
in  the  exploration  of  simpler  systems 
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